•'  Y » 


AEDC-TR-77-72 


JUL  2 2 3/7 


vhi 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  SEPARATED 
TURBULENT  FLOWS  ON  AXISYMMETRIC  AFTERBODIES 


NIELSEN  ENGINEERING  & RESEARCH,  INC. 
510  CLYDE  AVENUE 
MOUNTAIN  VIEW,  CALIFORNIA  94043 


June  1977 


Final  Report  for  Period  15  March  1976  to  15  April  1977 


Approved  for  public  release;  distribution  unlimited. 


Prepared  for 


°Pen-y  cr 


DIRECTORATE  OF  TECHNOLOGY 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 

AIR  FORCE  SYSTEMS  COMMAND 

ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


NOTICES 


When  U.  S.  Government  drawings  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a definitely  related  Government  procurement 
operation,  the  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

Qualified'  users  may  obtain  copies  of  this  report  from  the  Defense 
Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be 
considered  in  any  sense  as  an  endorsement  of  the  product  by  the  United 
States  Air  Force  or  the  Government. 


This  final  report  was  submitted  by  Nielsen  Engineering  & Research,  Inc., 
510  Clyde  Ave.,  Mountain  View,  California  94043,  under  contract 
F40600-76-C-0009,  with  the  Arnold  Engineering  Development  Center,  Arnold 
Air  Force  Station,  Tennessee  37389.  Mr.  Elton  Thompson  was  the  AEDC 
Project  Engineer. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 


APPROVAL  STATEMENT 

r 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 


ELTON  R.  THOMPSON 


Research  and  Development 
Division 

Directorate  of  Technology 


0 


ROBERT  O.  DIETZ 
Director  of  Technology 


UNCLASSIFIED 


| REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

IE 

2.  GOVT  ACCESSION  NO. 

1.  RECIPIENT'S  CATALOG  NUMBER 

4 T,Tte  (and  Sub“"m>  COMPUTER  PROGRAM  FOR  CALCU- 
LATION OF  SEPARATED  TURBULENT  FLOWS  ON 

S.  TYPE  OF  REPORT  A PERIOD  COVERED 

Final  Report  - March  15,  1976 
to  April  15.  1977 

AA1D  X IVUVlfj  1IUV  nr  X IiUDUL/lIj  Q 

8-  PERFORMING  ORG.  REPORT  NUMBER 

NEAR  TR  131 

7.  author^*; 

Gary  D*  Kuhn 

• . CONTRACT  OR  GRANT  NUMBERl'*)  . 

F40600-76-C-0009  \ 

• . PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Nielsen  Engineering  & Research,  Inc. 

510  Clyde  Avenue 

Mountain  View,  California  94043 

to.  program  element,  project,  task 
AREA  A WORK  UNIT  NUMBERS 

Program  Element  65807F 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Arnold  Engineering  Development  Center  (DYFS) 

12.  REPORT  DATE 

June  1977 

Air  Force  Systems  Command 

Arnold  Air  Force  Station.  Tennessee  37389 

11.  NUMBER  OF  PAGES 

92 

l«.  MONITORING  AGENCY  NAME  A AODRESSfff  dtlfmrmnt  Inm  Controlling  Olll cm) 

IS.  SECURITY  CLASS,  (ol  Ihlm  import) 

UNCLASSIFIED 

mwm  n^3|iM!H.>rTnr-Lr.viri'*™i  I 

IS.  DISTRIBUTION  STATEMENT  (ol  Utlg  Roporl) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  Ihm  (Attract  mind  In  Block  30,  If  dllftront  Dam  Rmport) 

*6.  supplementary  notes 

Available  in  DDC 

19.  KEY  WOROS  fConrJnu*  on  r*v*ram  aid*  11  nacaaaary  and  Idantlfy  by  block  number) 

compressible  flow  computer  programs 

finite-difference  methods  integral  methods 
interaction  methods  turbulent  boundary  layers 

separated  flow  transonic  flow 

inviscld  flow 

20.  ABSTRACT  (Contlnuo  on  rovormo  oldo  fi  nmcmmtory  and  Idontlfy  by  block  mtmbor) 

A computer  code  for  a turbulent  boundary-layer,  mviscld  interaction  method  for  axi- 
symmetric  configurations  of  the  type  used  for  isolated  nozzle  afterbody  models  Is  pre- 
sented. The  method  is  applicable  to  flows  with  subsonic  free  streams,  including  slightly 
supercritical  flows.  The  method  consists  of  an  integral  boundary-layer  method  and  a 
finite-difference  inviscid-flow  method  which  are  coupled  iteratively  through  the  boundary- 
layer  displacement  thickness.  Both  attached  and  separated  boundary  layers  can  be 

DD  I JAN  71  1473  EDITION  OF  I NOV  AS  1$  OBSOLETE 


UNCLASSIFIED 


UNCLASSIFIED 


20.  ABSTRACT  (Continued) 

calculated.  An  option  is  provided  for  calculating  two-dimensional  boundary  layers.  The 
procedure  for  separated  flows  is  to  specify  the  displacement  thickness  of  the  boundary 
layer  and  calculate  the  free-stream  velocity  distribution  from  both  the  boundary-layer 
equations  and  the  inviscid-flow  equations.  The  separation  point  location  and  the  angle 
of  the  displacement  surface  are  found  by  an  iterative  procedure.  The  equations  pro- 
grammed are  presented  along  with  detailed  instructions  for  the  preparation  of  input  data, 
description  of  the  program  output  and  instructions  for  operation  of  the  program  on  an 
IBM  370  computer.  Sample  cases  are  provided  for  a complete  axi symmetric  interaction 
calculation  and  for  a two-dimensional  boundary-layer  calculation.  ‘ ' 

• -v  : j : . ■ 


UNCLASSIFIED 


PREFACE 


"•'The  results  reported  herein  were  developed  for  the 
Arnold  Engineering  Development  Center  by  Nielson  Engi- 
neering and  Research,  Inc.  under  Contract  F40600-76-C-0009 . 
The  author  of  this  report  was  Gary  D.  Kuhn.  The  Air  Force 
Project  Engineer  for  this  contract  was  E.  R.  Thompson, 
AEDC/DYR.  The  Program  Element  Number  was  65807F.  This 
report  covers  the  work  done  during  the  period  15  March  1976 
to  15  April  1977.  The  reproducibles  used  in  the  reproduc- 
tion of  this  report  were  provided  by  the  author. 


AE  DC-TR-77-72 


CONTENTS 

Page 

1.0  INTRODUCTION  . . . . 7 

2.0  DEVELOPMENT  OF  BOUNDARY -LAYER  METHOD  8 

2.1  ASSUMPTIONS  8 

2.2  BOUNDARY -LAYER  EQUATIONS  FOR  COMPRESSIBLE 

TURBULENT  FLOW 9 

2.3  TRANSFORMATION  OF  AXISYMMETRIC  BOUNDARY -LAYER 

EQUATIONS  TO  ALMOST  TWO-DIMENSIONAL  FORM  ...  11 

2.4  TRANSFORMATION  OF  THE  COMPRESSIBLE  BOUNDARY- 

LAYER  EQUATIONS 12 

2.5  DEVELOPMENT  OF  INTEGRAL  BOUNDARY-LAYER  METHOD  . 15 

2.5.1  Integral  Equations  ......  15 

2.5.2  Velocity  Profiles 16 

2.5.3  Eddy  Viscosity  16 

2.5.4  Transitional  Eddy  Viscosity 17 

2.5.5  Equations  Solved  18 

2.5.6  Method  of  Solution  of  Equations 20 

3.0  INV  ISC  ID -FLOW  MODEL.  . . . 20 

4.0  VISCOUS -INVISC ID  INTERACTION  METHOD 22 

4.1  ESTIMATION  OF  SEPARATION  POINT  LOCATION  ....  23 

4.2  ESTIMATION  OF  BOUNDARY  OF  SEPARATED  REGION.  . . 24 

;4.3  ITERATION  PROCEDURE  25 

5.0  COMPUTER  PROGRAM  ORGANIZATION 26 

6.0  INPUT  TO  THE  PROGRAMS 29 

6.1  TABULAR  FORM 29 

6.2  DICTIONARY  OF  INPUT  VARIABLES 30 

7.0  PROGRAM  OPTIONS.  39 

7.1  BOUNDARY -LAYER  OPTION 39 

7.2  INVISC  ID-FLOW  OPTION 40 

7.3  VISC ID-INVISCID  ITERATION  OPTION 41 

7.4  BOUNDARY-LAYER  INITIAL  CONDITIONS 42 


3 


AEDC-TR-77-72 


,■  ...  Page 

8.0  PROGRAM  OUTPUT  . . . ...,  _.T  43 

8.1  STANDARD  OUTPUT  . . . ."  . . 43 

8.2  SPECIAL  OUTPUT  MESSAGES ■ r 47 

9.0  PROGRAM  OPERATING  PROCEDURE *.  .-50 

\ • J, 

9.1  GENERAL  JOB  CONTROL  SEQUENCE 51 

9.2  JOB  CONTROL  EXAMPLES . '51 

9.2.1  Creating  Partitioned  Data  Sets  .....  52 

9.2.2  Starting  an  Iteration  Sequence  . . . 52 

9.2.3  Restarting  an  Iteration  Sequence  . . . •■.  53 

9.2.4  Executing  the  Boundary-Layer  Program 

Alone ■ 54 

9.2.5  Executing  the  Inviscid  Program  Alone  -j^-.  55 

10.0  NUMERICAL  EXAMPLES  '.  55 

10.1  AXISYMMETRIC  INTERACTION 56 

10.2  EXAMPLE  OF  RESTARTING  AN  INTERACTION 

CALCULATION 57 

10.3  TWO-DIMENSIONAL  BOUNDARY  LAYER  ........  57 

10.4  AXISYMMETRIC  INVISCID  FLOW  . . . ...  . • .'  :^5g 

REFERENCES 59 

TABLES 

I.  RELATION  BETWEEN  EXTERNAL  DATA  SETS  AND  INPUT/OUTPUT 

LOGICAL  UNIT  NUMBERS 61 

II.  INPUT  DATA  CARDS 62 

ILLUSTRATIONS 

Figure 

1.  Boundary-layer  coordinate  system . 64 

2.  Schematic  of  iteration  procedure 65 

3.  General  relationship  of  programs  and  data  files.  . . 66 


4 


AEDCTR-77-72 


• Figure  pa9e 

4.  Axisymmetric  body  for  sample  calculation  ...  . . . ■ 67 

5.  Input  data  for  interaction  calculation  on  body  of 

figure  4.  . . . . . . . - . 68 


6.  Selected  output  for  calculation  of  viscid-inviscid 


interaction  on  boattailed  body.  70 

7 . Input  data  for  restarting  calculation  of  viscid- 

'inviscid  interaction  after  4 iterations 82 

8.  Two-dimensional  configuration  for  boundary-layer  . 

calculation  83 

9.  input  data  for.  two-dimensional  boundary-layer 

calculation  ......  84 

10.  Output  from  two-dimensional  boundary-layer 

calculation  ■, . . 86 

11.  input  data  for  inviscid-flow  calculation  on  body  of 

figure  4 . . 88 

NOMENCLATURE - 90 


5 


AEDC-TR-77-72 


1.0  INTRODUCTION 

One  of  the  critical  areas  in  the  design  of  both  aircraft 
and  missiles  is  the  interaction  between  propulsive  jets  and  the 
external  flow  over  the  aerodynamic  shapes  from  which  they  issue. 
The  drag  of  the  afterbody  and  exhaust  nozzle  can  be  a signifi- 
cant fraction  of  the  total  vehicle  drag.  In  order  to  obtain 
knowledge  of  the  afterbody  and  nozzle  flow  fields  and  their 
interaction  as  early  as  possible  in  the  design  procedure,,  it  is. 
desirable  to  have  a computational  method  for  accounting  for  the 
viscous  flow  over  the  afterbody  including  separation  and  reat- 
tachment of  the  flow.  Such  a method  should  also  account  for  the 
interaction  between  external  inviscid  flow  over  the  afterbody 
and  the  viscous  afterbody  flow  (including  separated  boundary 
layers) . 

This  report  describes  a method  for  predicting  the  viscous 
flow  field  about  an  axisymmetric  body  at  zero  angle  of  attack. 
The  method  combines  a finite-difference,  inviscid-f low  method 
with  an  integral  boundary— layer  method.  The  major  difficulties 
in  extending  existing  boundary-layer  technology  to  flow  on 
axisymmetric  afterbodies  are  the  same  as  those  encountered  in 
two-dimensional  flows.  First,  a model  for  the  turbulent 
Reynolds  stresses  must  be  employed.  For  the  kind  of  bodies  of 
interest  herein,  an  eddy-viscosity  model  has  been  found  which 
adequately  accounts  for  the  Reynolds  stresses.  For  afterbody- 
nozzle  configurations,  the  effect  of  the  jet  plume,  boundary- 
layer,  inviscid-f low  interaction  is  also  of  major  importance. 
Accounting  for  this  interaction  is  a difficult  problem.  In 
the  work  described  herein,  the  jet  plume  is  simulated  by  a 
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solid  body.  Thus,  the  boundary-layer  calculation  involves 
separation  on  the  boattail  section  of  the  afterbody  with 
subsequent  reattachment  on  the  simulated  plume. 

The  remainder  of  this  report  describes  briefly  the  equations 
which  are  the  basis  of  the  calculative  method  and  the  particular 
technique  used  to  calculate  the  viscous-inviscid  interaction. 

Then  the  computer  programs  are  described,  followed  by  detailed 
instructions  for  the  use  of  the  programs.  The  composite  computer 
program  described  herein  contains  the  boundary-layer  program 
developed  for  this  work  and  an  inviscid-flow  program  which  is 
essentially  the  same  as  that  described  in  reference  1.  Both 
programs  are  incorporated  as  subprograms  into  a mainline  program 
which  performs  an  iteration  to  calculate  the  interaction  between 
the  boundary  layer  and  inviscid  flow  for  axisymmetric  bodies  of 
the  type  used  for  nozzle  boattails  followed  by  a solid  plume 
simulator.  The  programs  are  written  in  the  FORTRAN  programming 
language  for  use  on  an  IBM  370  computer.  The  boundary-layer 
program  retains  the  capability  to  calculate  two-dimensional  ■ 
flows  as  an  option.  Both  the  boundary— layer  program  and  .the 
inviscid-flow  program  pan  be  used  individually,  without 
iterating  if  so  desired. 


2.0  DEVELOPMENT  OF  BOUNDARY -LAYER  METHOD 

The  complete  derivation  of  the  governing  equations  for  the 
boundary  layer  has  been  presented  in  reference  2.  In  this 
report  the  derivation  is  summarized  and  modified  for  the 
application  of  interest  herein. 

2.1  ASSUMPTIONS 

The  analysis  is  based  on  the  following  assumptions: 

1.  The  governing  equations  are  those  for  a compressible 
turbulent  boundary  layer. 
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2.  The  air  behaves  as  an  ideal  gas. 

3. -'  The  molecular  viscosity,  p.,  is  proportional  to  the 

temperature. 

'4.  The  specific  heat  of  the  gas  is  constant. 

5.  The  wall  is  either  two  dimensional  or  axisymmetric,  but 

can  have  arbitrary  profile  in  the  direction  of  flow  as 

. O'  *:  i 

long  as  the  longitudinal  radius  of  curvature  of  the 
wall  is  large  compared  to  the  boundary  layer. 

6.  The  pressure  is  constant  normal  to  the  wall. 

7.  The  wall  temperature  is  constant. 


2.2  BOUNDARY-LAYER  EQUATIONS  FOR 
COMPRESSIBLE  TURBULENT  FLOW 

The  basic  notation  and  coordinate  scheme  are  shown  in  fig- 
ure 1;^  Note  that  the  same  symbols  are  used  for  the  physical 
coordinates  of  both  two-dimensional  and  axisymmetric  configura- 
tions.; r Thus,  r denotes  the  distance  of  a point  from  the  axis 
of  an- axisymmetric  configuration,  or  from  the  reference  plane 
of  a two-dimensional  configuration,  x is  the  distance  along 
the  axis,  or  reference  plane  measured  from  the  leading  edge, 
and  the  dimension  y is  measured  from  the  body  surface  normal 
to  the  axis. 

The  governing  equations  describing  the  steady  flow  of  a 

■■  ?:  , 

compressible  turbulent  boundary  layer  for  both  two-dimensional 
and  axisymmetric  configurations  are: 


Continuity 

& (rV)  + ^ [rkp(v  - r^u)]  = 0 


Momentum 


p“-£*p<v-V»  I) 


(1) 


(2) 
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Energy 


dS  . < 1 x c^S 

pu  + p (v  . ^u) 


N 

O 


rk  by 


h + ir  (0-D 


6S 

Sy 


+ \x 


Kl-^)  + 4'  p. 


u ou 

jHe  3y 


(3) 


where  k = 0 for  two-dimensional  flow  and  k = 1 for  axisym- 
metric  flow,  r'  is 
respect  to  x,  and 


metric  flow,  is  the  derivative  of  the  body  radius  with 


S = 


- 1 


(4) 


Equations  (1) , (2) , and  (3)  are  easily  applicable  to  lami- 

nar and  transitional  flow.  In  laminar  flow,  substitution  of 
(3  = PT  = 1 reduces  the  equations  to  those  for  a laminar  bound- 
ary layer.  Further,  suitable  variation  of  the  eddy  viscosity 
and  turbulent  Prandtl  number  makes  the  equations  applicable  to 
the  transition  region. 

The  boundary  conditions  for  this  system  of  equations  are: 


y = 0: 


y = u = ue(x) 

du/dy  = 0 

v = ve  (x) 
S = 0 

x = xQ:  u = uQ(y) 

s = SQ(y) 
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2.3  TRANSFORMATION  OF  AXISYMMETRIC  BOUNDARY -LAYER"'  . 

EQUATIONS  TO  ALMOST  TWO-DIMENS IONAL  FORM 

in  order  to  put  the  axisymmetric  equations  (k  = 1)  into  a 
more  convenient  form,  the  Probstein-Elliott  transformation 
(ref.  3),,  is  applied.  The  coordinates  of  the  axisymmetric  body 
are  shown  in  figure  1.  The  Probstein-Elliott  transformation  is 


i y,G  ..* 

t i .£  >/  < 

~ rrw(x)r , 

dX  - [ L J dX 

(5) 

d?  - ZiZAL  dy 

(6) 

where  r (x)  is  specified  by  the  body  shape  and  r(x,y)  is 
w 

given • by  « 

' / r(x,y)  - rw(xj  + y (7) 

- ! r ;-:oo : . 

The  transformed  continuity  equation  has  the  form 

i Cii?  \jS  ’ # ‘ . 

d(P"?'  + k(.pv).  = o (8) 

- eno;  ;•  ox  by 


where 


and 


bib  rL  , ^ L ■'&. 

pV  = - = — P <v  “ rwu)  + T bi  pu 

dx  rw  rw 


(9) 

(10) 


Applying  the  transformation  to  the  momentum  and  energy  equa- 
tions (2)  and  (3)  yields  the  transformed  equations 
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pS  ^ + pv  3§L  - [(1-  + kty)  (mA  ^ + n if"  5 ^ 

dx  by  by  L \ dy  e by/ 


where 


A . i _ jl  + -a. 

P PT  PT 


B - p - £ + t ' p 


and  t is  the  transverse  curvature  factor 


t = a 
2 

rw 


(12) 


(13) 

(14) 


(15) 


and,  from  equations  (6)  and  (7) , 


y = 


(16) 


For  flows  in  which  the  transverse  curvature  terms  are 
negligible,  letting  k = 0 in  equations  (11)  and  (12)  produces 
the  equations  of  a two-dimensional  boundary  layer.  The  trans- 
verse curvature  terms  may  be  negligible  for  an  axisymmetric 
flow  if  the  body  radius  is  large  compared  to  the  boundary-layer 
thickness.  The  Probstein-Elliott  transformation  is  thus  a 

first-order  correction  of  the  approximate  equations  for  the 
effect  of  transverse  curvature,  allowing  the  boundary-layer 

thickness  to  be  of  the  same  order  as  the  body  radius. 


2.4  TRANSFORMATION  OF  THE  COMPRESSIBLE 
BOUNDARY-LAYER  EQUATIONS 

In  order  to  simplify  further  the  equations,  the  Stewartson 
transformation  (ref.  4)  reduces  the  equations  to  a set  of  equa- 

i 

tions  for  an  incompressible  flow.  The  following  analysis  is 
presented  in  terms  of  the  Probstein-Elliott  coordinates,  x and 
y,  with  the  understanding  that  in  the  two-dimensional  case 
x = x and  y = y. 
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In  the  Stewartson  transformation  the  following  variables 

r 

are  introduced: 


X 


■/ 


p a 
^e  e 


p a 

e 

o o o 


dx 


-I 


P a 

Fe  .e.  dy 

Pe  ae  Pe 
o o 


(17) 


„ o ~ 

U = — — u 


(18) 


V = 


S A f peac 

^ Peoa' 


pe  ae 

-P-  dy  + — — 

Pe  y Peae  P 


o p ~ 

1 v 


(19) 


With  these  the  boundary  conditions  become: 


y = 0: 

V “I  0 

VO^'O  ( *. 

.*•<:?  ?> 


U = V = 0 
S = Sw 

U = Ue  = aeQMe 


dU/dY  = 0 
S = 0 


It  is  assumed  that  S and  the  eddy-viscosity  parameter,  (3, 
transform  directly;  that  is, 


S(X,Y)  = S(x,?) 


(20) 


P(X,Y)  = P(x,y) 


(21) 


it  is  easily  shown,  using  relations  (17)  through  (21)  and  the 


perfect  gas  assumption  along  with  the  relations 


= c -3 

^ T, 


(22) 
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and 


= o 

dy 


(23) 


that  the  boundary-layer  equations  in  the  Stewartson  plane  are: 
Continuity 

I + ' (24) 


Momentum 


ri  au 

u 5x  + v Sy 


- (S  + D Ue  lx  * Cle  & [(1  +kt5)|i  0?]  <25) 


Ene rgy 


„ 5s  as 

U §X  + V dY 


- cve  57  [<1+kt?>A  if] 


+ c -srirte  [ ( 1 + kty)  BU  f ] 


(26) 


The  coordinate  y is  not  transformed  in  the  transverse  curva- 
ture terms  because  integration  of  the  equations  across  the 
boundary  layer  is  anticipated  and  only  corresponding  values  are 
needed  in  those  terms.  The  Chapman-Rubes in  parameter,  C,  is 
assumed  to  be  at  most  a function  of  x.  In  the  work  presented 
herein,  C is  a constant  evaluated  at  the  wall  temperature; 
that  is, 


where  Sutherland's  law  is  used  to  evaluate  the  viscosity. 
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M- 


T3/2 

T + Ts 


(28) 


where  A and  T are  constants, 
s 

In  the  remainder  of  this  report,  the  solution  of  the 
energy  equation  (26)  will  be  approximated  by  the  Crocco  relation 


S = 


(29) 


Thus,  the  velocity  profiles  found  to  be  valid*  for  incom- 
pressible', two-dimensional  turbulent  boundary  layers  can  be  used 
by  simply  transforming  the  input  quantities  to  the  incompres- 
sible plane,  performing  the  calculation  for  an  equivalent  incom- 
pressible boundary  layer,  and  then  transforming  the  results  back 
to  the  compressible  plane,  and  for  an  axisymmetric  flow,  back  to 
the  axisymmetric  coordinates. 


2.5  DEVELOPMENT  OF  INTEGRAL 
BOUNDARY -LAYER  METHOD 

2:5.1  ~ Integral  Equations 

The  integral  method  used  herein  was  described  in  detail  in 

> i ! i ! * ‘ ■ 

reference  5.  Families  of  integral  equations  are  derived  by 
eliminating  V between  the  momentum  and  continuity  equations 

' ■ ' . . ri "T  t.  ' 

and  then  taking  weighted  integrals  of  the  resulting  equation 

' j '<  u J : _ ■ 

across  the  boundary  layer. 

(30) 


(31) 


rH  r 1 

I UUX-UY  I Uxdn-  (S  + l)Ue(Ue)x-vOUy)Y  f(Y) 


dY  = 0 


In  the  present  case,  the  functions 


f = Yn  ; n = 0,1 
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produce  the  momentum  and  moment  of  momentum  integral  equations, 
respectively. 

2.5.2  Velocity  Profiles 

The  Y dependence  of  the  integral  equations  is  eliminated 
by  substituting  an  appropriate  parametric  formulation  for  the 
velocity  profiles.  The  function  used  for  the  present  theory  is 
a modification  of  Coles 1 family  (ref.  6)  with  a laminar  sublayer 
added  and  the  wake  function  approximated  analytically. 


U = U^[2.5  ln(l  +Y+)  + 5.1  - (3.39Y+  + 5 . l)exp  (-0 . 37Y+)  ] 


+ ViI#(f)37 


(32) 


The  parameter  is  the  friction  velocity, 


(33) 


The  variable  Y+  is  defined  to  account  for  the  axisymmetry  of 
the  flow 


(34) 


The  other  parameters  in  equation  (32)  are  <5^,  the  boundary- 
layer  thickness  and  a wake  velocity.  The  exponential 

terms  and  the  additional  unit  in  the  logarithmic  term  provide 
a smooth  transition  from  the  turbulent  flow  to  the  wall  through 
a laminar  sublayer. 


2.5.3  Eddy  Viscosity 

The  eddy-viscosity  model  used  in  this  work  is  an  extension 
of  the  two-layer  model  used  by  Kuhn  (ref.  2)  including  an 
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intermittency  function  for  the  outer  layer  and  a modification 
of  the  outer  layer  for  adverse  pressure  gradients  and  separated 
flows.  In  the  inner  layer  of  attached  flows,  the  eddy-viscosity 
parameter,  (3,  is  represented  by  an  exponential  expression  based 
on  the  law  of  the  wall.  In  the  outer  layer  Clauser's  expression, 
modified  for  adverse  pressure  gradients,  is  used  along  with  an 
intermittency  function  giving 

3=  [0.013  + 0.0038  exp(-6*Px/15Tw) ]Ue6j/[l + 5.5{y/6)6]  (35) 

For  favorable  pressure  gradients,  the  exponential  term  in 
equation  (35)  is  taken  to  be  unity. 

For  separated  flows,  the  eddy  viscosity  across  the  entire 
layer  is  represented  by  a relation  based  on  the  velocity 
profile  above  the  U = 0 line. 

U 5 

(3  = 0.013[1  + 5.5(y/6)6]-1  -f-  J (l  - ~ ) dy  (36) 

yu=o 

2.5.4  Transitional  Eddy  Viscosity 

Transition  from  laminar  to  turbulent  flow  is  calculated 
by  letting  the  eddy  viscosity  change  from  a laminar  viscosity 
to  a fully  turbulent  value  over  a short  distance  according  to 
the  relation 

3t  = [1  - exp  ( -K  (x  - xfc)  2)  ] Ot  - 1)  + 1 (37) 

where 

p is  the  transitional  eddy  viscosity 

PT  is  the  turbulent  eddy  viscosity 

x is  the  location  of  the  beginning  of  transition 
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and 


K = 


0.0001 


where 


6^  is  the  boundary-layer  thickness  at  x.  . 
t *- 


2.5.5  Equations  Solved 


Substitution  of  equation  (32)  into  the  two  equations  pro- 
duced by  equations  (30)  and  (31)  produces  two  ordinary  differ- 
ential equations  for  the  variation  of  the  variables  U^,  U^, 

6.,  and  U with  x.  A third  equation  produced  by  evaluating 
i ’ e 

equation  (32)  at  Y = 6^  allows  the  elimination  of  from 

the  equations,  leaving  a set  of  two  equations 


*11  (Vx  + Al26i  + Al3(Ue>x  ■ -UxlDTl/Ue6i 

X 


*ai(DT>x  + A336i  +A33<Ue>x  = 

X 


h 

■ts 


puy  dY 


e i o 


(38) 


(39) 


Anticipating  the  development  of  a viscous-inviscid  interaction, 
the  velocity,  Ue,  is  considered  to  b,e  a dependent  variable.  The 
coefficients  A^j  are  functions  of  the  variables  U^,  6^,  and 
Ue.  The  usual  procedure  for  solving  equations  (38)  and  (39) 
for  attached  boundary  layers  is  to  prescribe  the  pressure  dis- 
tribution or  the  free-stream  velocity  distribution,  Ue. 

However,  if  separation  occurs,  the  pressure  distribution  cannot 
be  prescribed  arbitrarily  in  the  separated  region.  If  an 
adverse  pressure  gradient  is  prescribed  for  an  attached  bound- 
ary layer,  the  value  of  UT  can  approach  zero.  When  UT 
vanishes,  the  coefficients  Ai:l  and  &22  in  equations  (38) 
and  (39)  also  vanish,  producing  a singularity.  The  singularity 
can  be  removed  by  rearranging  the  equations  so  that  U,f  is  not 
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a dependent  variable.  One  method  of  accomplishing  this  is 
simply  to  rearrange  the  equations  so  that,  can  be  prescribed 

and  Ug  can  be  calculated  as  a dependent  variable  as  shown  in 
reference  7. 

Another  method  of  avoiding  the  singularity  at  = 0 is 
used  in  the  method  described  herein.  The  displacement  thickness 
6*,  is  expressed  in  terms  of  u^,  6^,  and  ue  and  the  result  is 
differentiated  with  respect  to  x,  producing  a third  equation. 


; A (U  ) + A S . + A ( U ) + u 5*  =•  0 

31  T x 32  1 . 33  e x ex 


(40) 


where  dthe  A . are  functions  of  u , 5/,  and  n and  A 

3 1 T * 1 * e 31 

vanishes  when  U_  = 0.  It  can  be  shown  that  A /A 

x ir  3i 

and  A£ 1/A3  are  both  finite  when  u = 0.  This  allows  the 
three  equations,  (38),  (39),  and  (40),  to  be  reduced  to  two 

ordinary  differential  equations  in  the  four  dependent  variables, 
UT,  6*,  6^,  and  Ue.  The  resulting  equations  are: 


A«  Vx  +(a31  A32-a12)  51x+(a31  A==_,0  <ue'x 


u u 

T T 

U 6. 
e l 


J . . I • 


(41) 


A t / A„ , \ 

U'.6*:-+  ( ~~  a3p  -a„p  ) 6 . 
A3i  te  x VA31  32  £V  lx 

r. u : 


kJ7a33-a23)(V-  = 


= -(■ 

* . us2  J 


bUydY  (42) 


e i o 


- *•  . ;;r.  cv  f.  f.  . . 

The  value  of  can  be  obtained  directly  from  known  values  of 

6*,  6 • y " and  Ue  by  solving  the  nonlinear  relation 


6*  - £(UT,6i,Ue) 


(43) 
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Prescribing  the  distribution  of  6*  in  equations  ( 41 ) 
and  (42)  and  obtaining  UT  from  equation  (43)  is  equivalent 
to  prescribing  UT  as  described  previously. 

2.5.6  Method  of  Solution  of  Equations 

The  equations  (38)  and  (39),  or  (41)  and  (42)  are  integrated 
numerically  using  a fourth-order  Adams-Moulton  predictor-correc- 
tor integration  with  a fourth-order  Runge-Kutta  scheme  used  to 
obtain  starting  values.  By  doubling  or  halving  the  integration 
step  size,  the  integration  is  capable  of  some  optimization  of 
the  integration.  Thus,  numerical  errors  are  kept  within  certain 
bounds  by  dividing  the  integration  step  by  2 whenever  the  error 
is  too  large  and  multiplying  the  interval  by  2 whenever  less 
accuracy  is  necessary  than  is  being  obtained.  The  step  size  is 
allowed  to  decrease  to  as  small  a value  as  10~6  before  the 
calculation  is  terminated. 

3.0  INVISCID-FLOW  MODEL 

The  inviscid-flow  method  employed  herein  is  the  finite- 
difference  solution  of  the  full  potential  equation  described 
in  references  1 and  8.  The  computer  program  used  is  that 
described  in  reference  1 with  minor,  modifications  to  accommo- 
date the  iteration  between  the  inviscid  and  boundary-layer 
flows.  The  detailed  derivation  of  the  flow  equations,  the  cal- 
culative  method,  and  the  computer  program  are  contained  in 
references  1 and  8 and  will  not  be  repeated  here.  A general 
description  of  the  program  is  summarized  here  from  reference  1. 

"One  of  the  important  considerations  when  trying  to  solve 
the  full  potential  equation  is  the  choice  of  a coordinate  system. 
For  complex  three-dimensional  shapes  cartesian  coordinates  may 
be  best;  however,  for  simpler  two-dimensional  or  axisymmetric 
shapes  the  use  of  a coordinate  transformation  such  that  the 
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body, lies  along  a coordinate  line  can  greatly  simplify  the 
application  of  the  exact  boundary  condition  at  the  body  surface. 
The  program  described  in  this  paper  uses  a body-normal  coordi- 
nate system  for  closed  bodies.  For  open  bodies  {i.e.,  bodies 
with  a sting  or  simulated  wake)  it  uses  a body-normal  system  on 
the  forebody  up  to  the  first  horizontal  tangent  and  a sheared 
cylindrical  coordinate  system  aft  of  that  point.  This  coordi- 
nate . system  is  suitable  for  closed  bodies  which  are  blunt  on 
both  ends  and  convex  and  smooth  over  the  entire  body  or  for 
open  bodies  which  are  blunt-nosed  and  convex  and  smooth  up  to 
the  first  horizontal  tangent.  It  is  possible  to  treat  pointed 
bodies  and  bodies  with  slope  discontinuities  but  the  coordinate 
system  is  not  well  suited  for  them  and  their  solution  may  not 
be  as . accurate  as  the  blunt-body  solutions." 

"A  stretching  is  applied  to  both  the  normal  and  tangential 
coordinates  such  that  the  infinite  physical  space  is  napped  to 
a finite  computational  space.  Thus,  the  boundary  condition  at 
infinity  can  be  applied  directly  and  there  is  no  need  for  an 
asymptotic  far-field  solution.  Details  about  the  stretching 
functions  are  given  in  [reference  1]  appendix  A." 

"The  general  method  of  solution  is  to  replace  the  governing 
second-order  partial  differential  equation  with  a system  of 
finite-difference  equations,  including  Jameson’s  rotated  differ- 
ence scheme  at  supersonic  points.  The  difference  equations  are 
solved  by  a column  relaxation  method." 

"The  boundary  condition  at  the  body  surface  is  applied 
through  the  use  of  dummy  points  inside  the  body.  Details  of 
this  computation  are  given  in  [reference  1]  appendix  B." 
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4.0  VISCOUS -INVISC ID  INTERACTION  METHOD 

Separation  is  the  result  of  an  adverse  pressure  gradient 
causing  reversal  of  the  low-energy,  low-momentum  fluid  near  the 
wall  in  the  boundary  layer.  The  adverse  pressure  gradient  is 
determined  by  the  inviscid  flow  field  and  is  the  net  result  of 
the  body  shape  and  the  displacement  effect  of  the  boundary  layer . 
Thus,  the  viscous  and  inviscid  flow  fields  are  interdependent. 
Accounting  for  this  interdependence  with  an  interaction  model 
is  made  difficult  by  the  fact  that  the  boundary-layer  equations 
are  parabolic  and  therefore  cannot  respond  to  disturbances  down- 
stream of  local  stations,  while  for  subsonic  and  transonic  flows 
the  inviscid  flow  is  elliptical  and  therefore  subject  to  influ- 
ence by  the  entire  flow  field.  In  the  method  described  herein, 
the  boundary-layer  and  inviscid-flow  methods  are  used  alternately 
in  an  iterative  scheme. 

Before  proceeding  with  the  description  of  the  iterative 
scheme  developed  in  the  work,  it  is  worthwhile  to  consider  two 
important  results.  The  work  of  Presz  {ref.  9)  has  shown  that 
the  displacement  effect  of  the  separated  boundary  layer  on  an 
afterbody  can  be  represented  by  a conical  surface  placed  between 
the  separation  and  reattachment  points.  Reubush  and  Putnam 
(ref.  10)  went  a step  further,  and  devised  a calculative  scheme 
in  which  the  conical  surface  was  treated  as  a solid  surface  and 
the  boundary  layer  was  calculated  on  the  modified  body.  No 
reversed  flow  was  calculated  and  the  separation  prediction 
method  used  in  that  work  was  based  only  on  the  inviscid  pressure 
distribution  calculated  for  the  real  body  alone  and  thus  could 
not  account  for  Reynolds  number  effects.  Nevertheless,  these 
methods  reveal  significant  information  about  the  nature  of  the 
boundary  layer  in  the  region  of  strong  interaction.  First,  the 
boundary  layer  grows  smoothly  as  it  encounters  the  adverse 
pressure  gradient,  and  second  the  displacement  surface  assumes 
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a condition  such  that  its  effect  after  separation  can  be 
approximated  by  a conical  surface. 


4.1  ESTIMATION  OF  SEPARATION  POINT  LOCATION 


The  first  step  in  calculating  the  viscous-inviscid  inter- 
action is  to  calculate  the  inviscid  flow  over  the  plain  body. 
The  resulting  distribution  of  the  velocity  at  the  boundary,  uQ, 
is  then  prescribed  for  the  boundary-layer  calculation  using 
equations  (38)  and  (39).  If  strong  adverse  pressure  gradients 
exist,  the  bound ary -layer  calculation  may  reach  a point  where 
the  skin-friction  coefficient  approaches  zero  and  the  numerical 
calculation  can  proceed  no  further  with  ue  prescribed.  This 
point  is  not  necessarily  the  true  location  of  separation,  how- 
ever, since  the  interaction  with  the  inviscid  flow  has  not  yet 
been  accounted  for. 

A better  first  approximation  to  the  location  of  separation 
xg  can  be  found  by  computing  the  shape  factor,  Htr. 


(44) 


For  the  velocity  profiles  given  by  equation  (32)  this  parameter 
has  a value  of  4.0  when  U =0.  However,  experimental  measure- 
ments indicate  that  separation  actually  occurs  when  Htr  is 
approximately  2.0.  This  suggests  that  the  present  boundary- 
layer  approximation  is  not  accurate  in  the  vicinity  of  the 
separation  point.  The  major  cause  of  the  inaccuracy  is  believed 
to  be  the  neglect  of  any  upstream  influence  of  the  separated 
region  inside  the  boundary  layer.  However,  the  effect  of  the 
separated  region  which  is  transmitted  upstream  in  the  boundary 
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layer  is  believed  to  be  small  compared  to  the  effect  transmitted 
through  the  inviscid  flow.  In  this  work,  the  first  approximation 
to  the  value  of  xs  is  the  location  at  which  Htr  = 1.5  for 
the  inviscid  velocity  calculated  for  the  plain  body.  , 

4.2  ESTIMATION  OP  BOUNDARY  OF  SEPARATED  REGION 

At  the  separation  point,  xs,  the  skin  friction  is  assumed 
to  have  become  zero  and  the  calculation  is  carried  on  into  the  N 
separated  flow  using  a prescribed  distribution  of  5*  with 
equations  (41)  to  (43) . The  result  of  that  calculation  is  a 
solution  for  the  free-stream  velocity  (the  "viscous  velocity") 
which  may  or  may  not  agree  with  the  "inviscid  velocity"  produced 
by  calculating  the  inviscid— flow  theory.  An  iterative  procedure 
is  used  to  find  the  particular  variation  of  6*  downstream  of 
xs  for  which  the  "viscous  velocity"  and  the  "inviscid  velocity" 
agree. 

A two-parameter  analytical  formulation  is  used  to  represent 
the  effective  displacement  surface  between  separation  and  a 
point  downstream  of  reattachment.  Thus,  the  effective  body 
shape  is  given  by 

r = r + 6*  for  0 < x < x^  and  x < x < °°  (45) 

w s,  P 

r = (rw  + 6*)s  + (x  - xs)tan  Q for  x < x < x (46) 

where  Xp  is  the  location  of  the  peak  inviscid  pressure  as 
given  by  the  previous  inviscid  solution  and  03  is  the  angle 
of  the  extrapolated  6*  surface  with  the  axis. 

A first  approximation  to  the  angle,  0S,  is  obtained  from 
an  expression  presented  in  reference  10. 

9S  = tan-1  14.4  - 4.89  Ms  (47) 
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Downstream  of  the  free-stream  velocity  is  again  prescribed 

on  the  boundary  layer  (eqs.  (38)  and  (39))  as  an  exponential 
fairing  between  the  value  at  x^  and  the  inviscid  velocity 
from  the  previous  calculation.  ■ 


4.3  ITERATION  PROCEDURE 


The  method  developed  in  this  work  consists  of  three  itera- 
tions. ^One  iteration  is  used  to  locate  the  separation  point, 
xs.  Another  iteration  is  used  to  determine  the  angle,.;  0S,  of 
a conical  displacement  surface.  The  third  iteration  calculates 
the  best  solution  for  specific  values  of  xs  and  0_.  The 
iteration  procedure  is  described  schematically  in  figure  2. 

The  iteration  procedure  consists  of' three  cycles.  in  the 
inner  cycle/  the  inviscid  flow  and  the  boundary  layer  are  cal- 
culated alternately  until  the  largest  change  in  the  6*  solu- 
tion between  iterations  becomes  smaller  than  a specified  percent 
At  each 'step  of  the  cycle,  the  boundary-layer' displacement 
thickness  is  used  to  calculate  an  augmented  body  shape  by  the 
relation 

rn  * rw  + + (1'a)6^-1  ■ (48) 

where  rn  is  the  effective  body  radius  at  iteration  n and 
a is  a damping  factor,  usually  equal  to  0.5. 

When  the  inner  cycle  has  been  terminated,  the  calculation 
is  complete  if  no  separation  occurred.  However,  if  separation 
is  present,  there  exist  two  solutions  for  the  free-stream 
velocity  as  described  previously.  The  next  step  of  the 
procedure  is  to  calculate  the  squared  deviation  between  the, 
two  solutions  downstream  of  x_  at  19  points. 
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(49) 


where  is  the  "viscous  velocity"  and  uSi  is  the  "inviscid 

velocity."  This  quantity  is  then  compared  with  the  value  from 
the  previous  0S  iteration  and  if  a minimum  has  not  been  found, 
6 is  adjusted  in  an  appropriate  manner,  the  boundary  layer  is 
recalculated  using  the  new  value  of  9S  and  the  calculation 
reenters  the  inner  cycle.  This  process  is  repeated  until  a 
minimum  of  the  squared  deviation,  s,  is  found.  The  present 
program  calculates  the  best  fit  for  changes  in  9S  within 
0.5°. 


When  the  0g  cycle  is  complete,  for  a given  value  of  xs, 
the  next  step  is  to  adjust  xs  in  an  appropriate  manner  and 
reenter  the  6S  cycle.  These  cycles  are  repeated  until  the 
smallest  value  of  s is  found,  establishing  0g  within  0.5° 
and  xs  within  a distance  equal  to  the  value  of  5*  at  xs. 

An  alternate  termination  device  is  incorporated  in  the 
computer  program  to  aid  in  keeping  the  cost  of  the  solution  at 
a minimum.  If  the  squared  deviation,  s,  is  smaller  than  0.0005 
at  any  step  of  the  iteration,  the  iteration  procedure  terminates. 


5.0  COMPUTER  PROGRAM  ORGANIZATION 


In  this  section,  the  general  organization  of  the  programs 
will  be  described.  Specific  information  on  data  required  for 
input  and  data  developed  for  output  will  be  described  in  sec- 
tions 6,  7,  and  8.  Sample  Job  Control  Card  decks  for  a typical 
IBM  370  installation  are  presented  in  section  9. 

The  overall  program  consists  of  a mainline  program 
(program  1)  and  two  main  subprograms  each  consisting  of 
several  subroutines.  The  first  main  subprogram  is  the  inviscid- 
flow  program  (program  2) . It  is  a modified  version  of  the 
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program  described  in  reference  1.  The  second  main  subprogram 
is  the  boundary-layer  program  (program  3) . It  is  based  on  the 
integral  theory  described  previously  and  retains  the  two- 
dimensional  case  as  an  option.  The  mainline  program  controls 
the  iteration  between  the  other  two  programs.  Either  of  the 
two  subprograms  may  be  used  separately,  without  iterating  by 
appropriate  choice  of  the  input  parameters. 

Both  the  boundary-layer  program  and  the  inviscid-flow 
program  require  some  punched  card  input  and  some  input  data  from 
disc  or  tape  data  files.  The  data  files  must  be  identified  by 
specific  Logical  Unit  numbers.  Each  program  in  turn  produces 
new  data  files  and  printed  output.  The  general  relationship  of 
the  programs  and  the  various  data  files  are  shown  in  figure  3. 
The  specific  Logical  Unit  numbers  required  for  input  and  output 
are  listed  in  Table  I.  Two  Logical  Unit  numbers  are  associated 
with  each  data  file  shown  in  figure  3.  One  unit  is  used  for 
input,  the  other  for  output. 

Program  2 requires  initially  data  from  cards  describing 
the  free-stream  conditions,  the  computational  mesh  and  the  body 
shape.  Alternately,  the  program  can  accept  the  input  body  shape 

t 

from  data  file  1.  It  can  also  accept  an  initial  solution  for 
the  perturbation  velocity  potential  from  another  data  file  (data 
file  2) . Program  2 produces  printed  output  lists  of  the  appro- 
priate flow  field  quantities,  quantities  describing  the  configu- 
rations and  the  computational  mesh  and  several  data  files.  Data 
file  2 is  rewritten  using  the  new  solution  for  the  potential.  A 
third  file  (data  file  3}  is  written  containing  the  distribution 
of  the  axial  velocity  component  for  use  by  the  boundary-layer 
program. 

Program  3 requires  initially  the  free-stream  conditions  and 
gas  constants  as  well  as  parameters  describing  the  shape  of  the 
surface  over  which  the  boundary  layer  is  flowing.  On  the  first 
iteration  of  a viscous-inviscid  interaction,  the  surface  shape 
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is  the  same  as  that  for  program  2.  On  subsequent  iterations, 
the  body  shape  for  the  boundary-layer  program  remains  the  same 
while  that  for  the  inviscid  flow  is  modified  by  the  addition  of 
the  boundary  layer.  Program  3 can  also  accept  data  from  data 
files  as  optional  input.  The  free-stream  velocity  distribution, 
u0,  can  be  input  from  data  file  3,  produced  by  program  2.  The 
distribution  of  the  body  shape  augmented  by  the  displacement 
thickness  can  be  input  from  data  file  4.  That  file  differs  from 
data  file  1 because  it  contains  the  raw  data  for  6*  + r.,  versus 
x as  calculated  by  the  boundary-layer  program  while  data  file  1 
contains  the  shape  adjusted  according  to  equation  (48)  and  inter- 
polated to  the  x stations  of  the  original  input  shape.  If 
program  3 is  being  used  separately  from  program  2 (i.e.,  without 
iterating) , two  additional  options  are  available  for  card  input. 
Either  the  free-stream  velocity,  u0,  can  be  input  as  mentioned 
previously,  or  the  displacement  thickness,  6*,  may  be  input. 

These  options  are  described  more  fully  in  sections  6.2  and  7.1. 

Program  3 produces  as  output  lists  of  the  boundary  layer 
and  flow  quantities  as  they  are  calculated  along  the  body.  In 
addition,  program  3 produces  an  updated  version  of  data  file  -4 
and  the  augmented  body  shape  (data  file  1)  required  by  program  2. 
Data  file  4 also  contains  a list  of  the  velocity  ratio  ue/ueQ 
corresponding  to  the  boundary  layer. 

Program  3 can  be  used  in  a two-dimensional  mode  if  desired. 
However,  calculation  of  a viscous-inviscid  interaction  can  only 
be  done  for  axisymmetric  cases  with  the  present  inviscid  program. 

All  card  input  pertaining  to  programs  2 and  3 is  input 
through  program  1.  Program  1 also  produces  a file  of  the  quan- 
tities needed  to  restart  the  calculation  if  the  calculation 
should  terminate  before  all  iterations  are  completed.  These 
data  are  stored  on  data  file  5.  Detailed  instructions  regarding 
restarting  are  presented  in  section  7.3. 
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6.0  INPUT  TO  THE  PROGRAMS 

The  data  required  by  the  programs  generally  fall  into 
three  categories:  (1)  geometrical  data;  (2)  flow  field  data; 

and  (3)  control  parameters.  The  control  parameters  are  indices 
for  specifying  options  and  iteration  counters.  It  will  be  noted 
by  comparison  with  reference  1 that  a number  of  input  quantities 
required  for  the  inviscid-flow  program  have  been  eliminated  in 
the  present  version.  This  has  been  done  by  incorporating  the 
calculations  required  to  obtain  some  of  the  quantities  into  the 
present  code  or  by  simply  defining  fixed  values  which  have  been 
found  to  be  successful.  Specifically,  a value  of  1.4  is  used 
for  the  initial  value  of  the  subsonic  relaxation  factor,  a value 
of  0.1  is  used  for  the  initial  value  of  the  supersonic  relaxation 
factor  and  a value  of  1.3  is  used  for  the  exponent  in  the  normal 
coordinate  stretching  function.  • Also,  it  is  assumed  that  the 
computational  grid  has  equal  step  sizes  in  both  coordinate 
directions  at  the  nose  of  the  body. 

The  general  requirement  of  the  input  data  is  that  the  tabu- 
lar lists  of  the  various  distributions  required  represent  smooth 
curves.-  This  is  especially  true  of  the  list  of  body  shape  co- 
ordinates. The  inviscid  program  uses  cubic  splines  to  fit  the 
input  coordinates,  so  those  coordinates  must  accurately  represent 
a smooth  curve  with  continuous  second  derivatives. 

6.1  TABULAR  FORM 

The  input  data  required  for  calculating  transonic  boundary- 
layer,  inviscid-flow  interactions  consist  of  several  punched 
cards  containing  parameters  describing  the  free-stream  flow 
conditions,  the  computational  mesh  for  the  inviscid  calculation, 
initial  values  for  the  boundary-layer  calculation,  and  certain 
options  that  are  available  in  the  programs.  A dictionary  of 
the  input  data  is  presented  in  the  next  section.  Table  II  shows 
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the  input  variables  as  they  are  to  be  punched  on  the  data  cards. 
More  detailed  explanation  of  the  requirements  for  the  inviscid- 
flow  program  are  presented  in  reference  1 and  are  not  repeated 
herein. 

6.2  DICTIONARY  OF  INPUT  VARIABLES 

The  variables  required  for  input  on  punched  cards  are 
defined  in  this  section  in  the  order  in  which  they  are  required. 
Additional  details  on  the  format  of  the  punched  data  are  given 
in  Table  II.  The  first  three  cards  of  any  input  data  deck  con- 
tain a description  of  the  case  being  calculated.  Any  or  all  of 
these  three  cards  may  be  blank,  but  all  three  are  required.  The 
remaining  variables  in  Table  II  are  as  follows: 

LPROG  Integer  indicating  level  of  calculation. 

= -1  Inviscid  flow  only. 

= 0 Inviscid-viscous  interaction. 

= 1 Boundary  layer  only. 

NRSTRT  Integer  indicating  whether  calculation  is  being  restarted 
to  continue  a previous  calculation.  Only  needed  if 
LPROG  =0. 

- 0 Start  from  zero.  Input  all  quantities  on  cards 
or  data  files  as  required. 

> 0 Restart.  Input  data  file  5 (Logical  Unit  11) 
containing  data  from  previous  iteration  plus 
all  other  input  data  files. 

NPRINT  Integer  indicating  quantity  of  output  to  be  printed 
(see  section  8.1). 

= 0 Minimum  output. 

= 1 All  output. 

NlMAX  Integer  convergence  criterion  for  inner  iteration  cycle. 

The  cycle  is  considered  to  be  converged  when  the  maxi- 
mum change  in  6*  between  iterations  is  less  than 
NlMAX  percent.  A value  of  7 is  good  for  most 
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NlMAX 

(cone. 

Nl 

N2 

N3 

IBL 

MIT2 


GAMMA 

AMINF 

IXY 


XO,YO 


, engineering  purposes'.  More  accuracy  is  obtained  with 
smaller  values  at  increased  cost.  However,  a maximum 
of  11  iterations  will  be  performed  in  the  inner  cycle 
in  any  case. 

Integer  iteration  counter  for  inner  viscid-inviscid 
iteration. 

Integer  iteration  counter  for  6Q  iteration  presently 
limited  to  a maximum  of  20. 

Integer  iteration  counter  for  xs  iteration  presently 
limited  to  a maximum  of  20. 

Integer  indicating  how  interaction  calculations  are  to 
begin. 

= 3 Start  with  inviscid  flow. 

= 0 Start  with  boundary  layer. 

Integer  number  of  iterations  to  be  executed  by  the 
inviscid  program  in  solving  its  relaxation  equations 
for  the  first  series  of  calculations  in  the  inner 
.iteration  cycle.  After  the  inner  iteration  cycle  has 
been  converged  for  the  first  time,  the  relaxation 
iteration  limit  is  set  equal  to  the  value  of  MIT  to 
be  input  subsequently;  A typical  value  for  MIT2 
is  20  for  interaction  calculations.  A larger  value 
may  be  better  for  some  transonic  cases. 

Ratio  of  specific  heats. 

Free-stream  Mach  number. 

Integer  number  of  values  of  coordinate  pairs,  XO,YO,  to 
be  input  for  inviscid  body  shape.  If  IXY  = 0,  the 
required  shape  must  be  input  from  data  file  1 (Logical 
Unit  14) . Maximum  value  is  100. 

Axial  and  radial  coordinates  of  body  shape  for  inviscid- 
flow  calculation,  2 per  card. 
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The  next  series  of  variables,  items  6 and  7 in  Table  II,  are  for 

the  inviscid-f low  program.  Detailed  information  on  how  these 

data  are  to  be  obtained  is  contained  in  reference  1. 

IMAX  Number  of  grid  lines  in  the  tangential  direction; 

1=1  is  the  forward  stagnation  line,  I = IMAX  is  the 
rear  stagnation  line  for  closed  bodies  and  downstream 
infinity  for  open  bodies.  For  each  grid  refinement 
IMAX  is  increased  such  that  IMAXNEW  = 2(IMAX0LD)  - 1. 
The  present  limit  on  IMAX  is  81.  Instructions  for 
changing  this  limit  appear  as  comments  in  the  program 
listing  (subroutine  ONEO) . 

JMAX  Number  of  grid  lines  in  the  normal  direction;  J = 1 

corresponds  to  an  infinite  distance  from  the  body  and 
J = JMAX  is  on  the  body.  The  same  formula  and  limit 
that  apply  to  IMAX  also  apply  to  JMAX. 

MIT  Maximum  number  of  iterations  (complete  relaxation  cycles) 
allowed  for  the  inviscid-f low  program  after  complete 
inner  cycle  or  when  using  the  inviscid-f low  program 
alone.  A typical  value  for  interaction  calculations 
is  same  as  MIT2.  The  values  of  MIT  and  MIT2  are 
combined  for  the  final  iteration  of  an  interaction 
calculation  so  that  a total  of  MIT+MIT2  relaxation 
cycles  are  performed. 

MHALF  Number  of  grid  refinements  to  be  done.  For  interaction 
calculations  a value  of  zero  should  be .used  with  a 
grid  fine  enough  for  adequate  resolution. 

KLOSE  Body  type. 

= 0 Open  body  (i.e.,  one  with  a sting  or  wake) . 

= 1 Closed  body. 

NPLOT  Plot  trigger;  NPLOT  = 1 causes  write  on  disc  for  input 
to  plot  routines  and  calling  of  plot  routines.,  (The 
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plot  routines  are  not  included  in  the  present  version 
so  NPLOT  = 0 should  be  put  in.) 

LREADP  Integer  indicating  whether  initial  estimate  of  potential 
distribution  is  to  be  input  from  data  file  2 (Logical 
Unit  13) . 

= 0 - no 
e 1 - yes 

DETAO  Step  size  of  the  normal  coordinate  at  the  body. 

XIXM  Value  of  the  computational  coordinate,  X,  at  the 

matching  point  of  the  two  stretching  functions  used  in 
the  finite-difference  scheme  (see  ref.  1) , for  open 
bodies  only.  Since  X varies  from  zero  to  one,  XIXM 
is  the  fraction  ^of  the  total  number  of  grid  points 
which  will  be  in  the  first  stretching  region  (ahead  of 
xm) . Usual  value  is  about  0.85.  The  exact  value  used 
should  be  chosen  so  that  one  grid  point,  usually  the 
one  at  x^  corresponds  to  the  boattail-sting  junction. 

XM  Axial  location,  xm,  (in  physical  coordinates)  of  the 

junction  (or  matching  point)  between  the  two  tangential 
stretching  functions,  for  open  bodies  only,  see  refer- 
ence 1.  Must  be  less  than  XO(IXY).  This  parameter  is 
used  to  concentrate  computational  mesh  points  in  a 
certain  region.  The  usual  approach  for  interaction 
calculations  is  to  let  xm  be  equal  to  the  length  of 
the  body  to  the  beginning  of  the  sting,  XBT. 

XIMl  Approximate  location  of  the  last  finite  value  of  the 
tangential  coordinate  for  open  bodies;  usually  about 
2 or  3 maximum  diameters  larger  than  the  length  of 
the  body,  XBT. 

XBT  Length  of  the  body  to  the  beginning  of  the  sting. 


The  remaining  input  variables  are  related  to  the  boundary-layer 
program  alone  or  to  the  viscous-inviscid  interaction  method. 
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XZNEW  Axial  location  at  which  boundary-layer  calculations 

will  begin  after  four  iterations  have  been  calculated. 

The  usual  procedure  is  to  start  the  boundary-layer 
calculations  close  to  the  nose  of  a long  body  at  XZ 
(see  item  15  in  Table  II)  and  then  after  four  iterations 
move  the  starting  point  to  XZNEW  (XZNEW  > XZ)  . In 
• subsequent  iterations,  the  boundary  layer  does  not 

change  for  X < XZNEW.  For  long  slender  bodies  with 
boattails,  XZNEW  can  usually  be  near  the  beginning  of 
the  boattail. 

DTHET  Angular  deviation  of  extraplated  6*  surface  from  the 

programmed  value  given  by  equation  (47).  Usual  starting 
value  is  0.0°. 

LSEP  Logical  variable  indicating  whether  the  location  of 
separation  (XSEP)  is  known  (.TRUE.)  or  not  (FALSE). 

XSEP  Axial  location  of  separation.  Put  in  only  if  LSEP  = TRUE. 

IOPT  integer  indicating  the  mode  of  the  calculation. 

=1  ug  is  to  be  input. 

= 2 6*  is  to  be  input. 

Put  in  a value  of  1 for  starting  an  interaction 
calculation. 

K Integer  indicating  whether  flow  is  axisymmetric. 

= 0 Two  dimensional. 

= 1 Axisymmetric. 

LVARl  Integer  indicator  for  method  of  input  of  ue  when 
IOPT  =1. 

= 0 Input  ue  (dimensional)  on  cards. 

= 1 Input  u_/u^  on  cards. 

e cq 

= 2 Input  u_/uQ  from  data  file  3 on  Logical 

Unit  12?  S° 
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LSHAPE 


Lie 


LDSTAR 


IUNIT 


LSHPBL 


Integer  indicating  option  for  calculating  all  initial 
conditions  (I.C.)  except  ue  (see  section  7.4). 

= 0 Input  initial  values  per  LIC. 

= 1 calculate  I.C.  for  flat  plate. 

*=  2 Calculate  I.C.  for  cylinder. 

= 3 Calculate  I.C.  for  cone. 

Integer  indicating  initial  condition  options  for 
IOPT  = 1 and  LSHAPE  = 0. 

= 1 Put  in  CPC1  and  DELTAl. 

- 2 Put  in  CFCl  and  DELSTl. 

Integer  indicating  whether  a file  of  6*  + r w is  to 
be  input . 

= 0 No  input. 

= 1 File  of  6*  + rv  versus  x is  required  on 
Logical  Unit  15  (data  file  4) . 

Integer  indicating  which  value  of  the  gas  constant, 

RGAS,  and  the  constants  in  Sutherlands  temperature- 
viscosity  relation  (eq.  (28))  are  to  be  used.  The 
choice  depends  on  whether  air  is  the  gas  being  calcu- 
lated and  the  units  of  the  input  quantities. 

= 1 Input  units  must  be  pounds,  feet,  seconds,  and  °R. 

= 2 Input  units  must  be  pounds,  inches,  seconds, 

and  °R. 

= 3 Input  units  must  be  newtons,  meters,  seconds, 
and  °K. 

For  another  gas,  or  other  units,  put  in  anything  for 
IUNIT  and  put  in  nonzero  values  of  VISC,  RGAS,  and 
SCON. 

Integer  indicating  whether  body  shape  is  to  be  input 
for  boundary-layer  calculations. 

= 0 XRL  and  RL  are  assumed  to  be  the  same  as  XO  and 
YO.  This  is  usually  the  case  when  starting  an 
interaction  calculation. 

= 1 XRL  and  RL  will  be  required. 
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note: 

NVAR 

XVAR, 

VAR 


EX 

PT 

TT 

TWONTT 

vise 


RGAS 


The  next  three  variables,  NVAR, XVAR, VAR,  are  only 

required  on  cards  if  LPROG  «=  1 and  LVARl  = 0 or  1." 

Integer  indicating  the  number  of  values  to  be  input  for 
the  prescribed  variable  (ue  or  6*).  Maximum  value 
is  100. 

Axial  location  and  value  of  prescribed  variable  as 
follows: 

IOPT  o 1 and  LVARl  = 0,  VAR  = ue  / 

I0PT  = 1 and  LVARl  = 1,  VAR  = Ue/ue 
IOPT  = 2,  VAR  =6* 

Reference  length.  Needed  if  input  data  lengths  are 
nondimensionalized.  If  lengths  are  dimensional,  put 
in  EL  = 1.0. 

Total  pressure,  pfc  (lb/ft2),  (lb/in2) , or  (newton/m2). 

Total  temperature,  Tt  (°R)  or  (°K)  . 

Ratio  of  body  surface  temperature  to  total  temperature, 

W ■ - 

Constant  X in  Sutherlands  formula  for  viscosity, 
equation  (28)  . if  one  of  the  programmed  values  is 
acceptable,  put  in  a value  of  0.0.  The  value  used 
will  then  be  determined  by  the  value : of  IUNIT  on  item 
number  11  as  follows: 

IUNIT  VISC 

1 2.27(10-8)  lb  sec/ft2  (OR)1/2 

2 1.5764 (10“10)  lb  sec/in  (Or)1/2 

3 1.4582 (10“6)  Newton  sec/m  (OR)1/2 

Gas  constant.  If  one  of  the  programmed  values  is 
acceptable,  put  in  a value  of  0.0.  The  value  used 
will  be  determined  by  the  value  of  IUNIT  as  follows: 
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RGAS  I UNIT  RGAS 

(cone.)  1 1716.0  ft2/sec2  °R 

2 247104.0  in2/sec2  °R 

3 286.96  m2/sec2  °K 

SCON  Constant  Ts,  in  Sutherlands  viscosity  law,  equation  (28) . 

If  one  of  the  programmed  values  is  acceptable,  put  in  a 
value  of  0.0.  The  value  used  will  be  determined  by  the 
value  of  IUNIT  as  follows: 

IUNIT  SCON 

1 198.6  °R 

2 198.6  °R 

3 110.333  °K 

DFACT  Relaxation  factor  for  adding  6*  to  body,  the  factor 
a in  equation  (48) . Usual  value  is  0.5.  If  a value 
of  0.0  is  put  in,  a value  of  0.5  will  be  used. 

XZ  Axial  location  of  beginning  of  boundary-layer 

calculation. 

RLEN  Axial  location  of  end  of  boundary-layer  calculation. 

Usually  at  least  one  maximum  diameter  larger  than  XBT. 

XT  Axial  location  of  transition  from  laminar  to  turbulent 

boundary-layer  flow. 

DXZ  Integration  step  size  at  the  beginning  of  the  boundary- 

layer  calculation.  The  usual  procedure  is  to  use  a 
small  value  for  a new  case  until  experience  indicates 
whether  a larger  value  would  work.  If  too  large  an 
initial  step  is  attempted,  the  first  step  of  the 
integration  procedure  may  calculate  a negative  value 
of  6i  and  the  calculation  will  terminate.  A typical 
value  would  be  (10“4) *XBT. 

DXP  Axial  interval  at  which  velocity  profiles  are  to  be 

printed.  If  a value  of  0.0  is  put  in,  no  profiles  are 
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DXP 

(cone. ) 
DXOUT 


HLIM 


CPC1 

DELTA1 

DELSTl 

UEl 

DUEDX 


printed.  In  any  case,  profiles  will  not  be  printed 
more  often  than  the  output  step  specified  by  DXOUT. 

Axial  interval  at  which  output  of  the  local  values  of 
UT,  6,  <5*,  Cj,  Uq/uq^,  Cp,  3-nd  are 

to  be  printed.  The  output  may  be  at  irregular  intervals 
since  the  output  location  is  determined  as  the  point 
where  the  local  value  of  X first  becomes  greater  than 
XOUT.  The  new  value  of  XOUT  is  then  calculated  as 
XOUT  = X + DXOUT.  This  value  should  be  selected  to 
produce  less  than  200  printed  output  steps  for  the 
entire  length  of  the  body. 

Limit  value  of  Htr  to  indicate  separation.  Input  of 
a blank  or  a value  , of  0.0  will  cause  a value  of  1.5  to 
be  used  as  discussed  in  section  4.1.  if  the  interaction 
iteration  has  difficulty  converging  due  to  a very  small 
separated  region,  it  may  be  necessary  to  increase  this 
value  slightly  in  order  to  calculate  a fully  attached 
boundary  layer.  A value  greater  than  4.0  will  allow 
the  calculation  to  proceed  to  a singularity  if  one 
should  occur.  The  calculation  wilJ  terminate  at  that 
point. 

Value  of  skin-friction  coefficient  at  initial 
boundary-layer  station  (see  section  7.4). 

Value  of  boundary-layer  thickness,  6 (compressible), 
at  initial  bound ary- layer  station  (see  section  7.4). 

Value  of  boundary-layer  displacement  thickness,  6*,  at 
initial  bound ary -layer  station  (see  section  7.4). 

Value  of  free-stream  velocity  at  initial  boundary-layer 
station. 

Value  of  free-stream  velocity  gradient  at  initial 
boundary-layer  station. 
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NR  Integer  number  of  values  of  XRP  and  RL  to  be  input  for 

body  shape.  If  NSHPBL  = 0,  this  is  assumed  to  be  the 
same  as  IXY.  Maximum  value  is  100. 

XRP,RL  Axial  and  radial  coordinates  of  body  shape  for  boundary- 
layer  calculation.  If  LSHPBL  = 0,  these  are  assumed 
to  be  the  same  as  XO  and  YO,  respectively.  For  two- 
dimensional  configurations,  these  represent  the  x and  y 
coordinates  of  a surface  measured  from  a reference  plane. 

7.0  PROGRAM  OPTIONS 

Several  optional  modes  of  calculation  are  available  through 
the  input  parameters.  A description  of  the  options  and  the 
corresponding  values  of  the  pertinent  parameters  follows. 

7 . 1 BOUNDARY -LAYER  OPTION 

To  use  only  the  boundary-layer  program,  put  in  the  three 
card  description,  then  the  first  value  on  the  fourth  input  data 
card  should  be: 

LPR0G=1  (see  Table  II) 

The  remaining  variables  shown  on  item  number  2 in  Table  II 
are  not  required.  The  remaining  cards  would  be  those  corres- 
ponding to  item  3 and  items  11  to  19  as  described  in  Table  II. 
With  this  option,  the  user  has  the  choice  of  specifying  either 
the  free-stream  velocity,  ue,  or  the  displacement  thickness, 

<5*,  through  the  variables  NVAR,  XVAR,  and  VAR  on  items  12  and  13. 
The  boundary-layer  calculation  can  be  restarted  at  any  station 
by  inputting  the  values  of  the  variables  at  that  station  as 
listed  in  the  output.  The  calculated  list  of  6*  + rw  and 
ue/ue0  will  be  written  on  Logical  Unit  10  (data  file  4)  when 
the  calculation  terminates. 
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Another  method  is  also  available  for  calculating  the  bound- 
ary layer  alone.  The  boundary-layer  step  of  a viscid-inviscid 
iteration  can  be  executed  separately.  The  appropriate  values 
on  the  fourth  card  would  be: 

LPROG  = 0 
NRSTRT  = 0 or  1 
NPRINT  o 0 or  1 
NlMAX  = 100 
Nl  =1 
N2  = 21 
N3  = 21 
IBL  = 0 

MIT2  (not  needed) 

This  option  requires  input  of  all  quantities  as  though  the 
iterative  sequence  were  to  be  completed,  with  the  values  just 
described,  only  the  boundary  layer  will  be  calculated  and  then 
the  run  will  be  terminated.  If  it  is  desired  to  continue  the 
iteration,  simply  make  NlMAX  equal  the  desired  convergence  per- 
centage, make  N2  and  N3  less  than  21,  and  put  in  an  appropriate 
value  for  MIT2.  The  free-stream  velocity  distribution  must  be 
provided  on  Logical  Unit  12  for  this  case.  All  other  optional 
inputs  are  the  user's  choice. 

7.2  INVISCID-FLOW  OPTION 

To  use  only  the  inviscid-f low  program,  put  in  the  three 
card  description,  then  put  in  LPROG  = -1  on  the  fourth  card 
(item  number  2 in  Table  II) . Of  the  remaining  values  on  that 
card,  only  NPRINT  is  required.  The  value  of  NPRINT  is  the  user's 
choice.  After  the  first  four  cards  only  the  data  for  items  3 
to  7 as  described  in  Table  II  and  section  6 are  required  for 
this  option. 
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7.3  ' VISCID-INVISCID  ITERATION  OPTION 

To  use  both  the  boundary-layer  and  the  ihviscid— flow  pro- 
grams iteratively,  put  in  LPROG  = 0 and  all  other  quantities 
as  appropriate.  Such  iterations  can  be  started  with  only  the 
body  shape  and  free-stream  flow  quantities  known  and  may  be 
restarted  to  continue  a prematurely  terminated  iteration. 
Several  options  are  available  to  the  user  for  restarting  an 
unfinished  iteration.  See  section  10.2  for  an  example  of 
restarting.  The  simplest  option  is  to  put  in  NRSTRT  =1  as 
the  second  value  on  the  fourth  input  data  card  and  to  provide 
the  required  input  data  files  on  Logical  Units  11,  12,  13,  14, ; 
and  15  (see  fig.  3 and  Table  I) . The  required  values  of  the 
remaining  variables  on  card  number  4 are  printed  periodically 
in  the  output.  See  section  8.1,  step  20  of  the  output  list. 

The  particular  set  of  values  needed  for  restart  is  the  last  set 
printed  before  termination.  The  only  other  data  required  for 
restarting  are  the  three-card  description  of  the  case.  The 
calculation  then  proceeds  from  where  the  previous  iteration 
stopped.  Another  method  of  restarting  would  be  to  omit  the 
restart  file  and  put  in  NRSTRT  = 0.  The  user  can  then  vary 
any  of  the  other  input  quantities,  using  the  data  files  or 
punched  cards  as  desired;  Note  that  the  calculation  terminates 
when  N2  or  N3  reach  a value  of  21.  The  value  of  Nl  increases 
continuously  throughout  the  calculation  while  the  value  of  N2 
is  reset  to  1 each  time  the  6S  cycle  is  converged.  If  the 
xe  iteration  cycle  has  not  converged  when  N3  reaches,  the  value 
of  21,  that  value  must  be  decreased  in  order  to  continue. 
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7.4  BOUNDARY -LAY ER  INITIAL  CONDITIONS 


Initial  values  of  boundary-layer  quantities  can  be  obtained 
in  several  ways.  The  user  can  obtain  values  of  the  skin-friction 
coefficient,  Cf,  and  either  the  boundary-layer  thickness,  6,  or 
the  displacement  thickness,  6*.  These  are  shown  on  item  16  in 
Table  II.  The  appropriate  values  LSHAPE  = 0 and  LIC  = l^or  2 
are  then  punched  on  the  card  corresponding  to  item  11.  For  the 
case  when  no  other  source  of  this  information  is  available, 
formulas  have  been  included  in  the  program  based  on  the  Blasius 
solution  for  laminar  boundary  layers  and  based  on  one-seventh 
power  law  velocity  profiles  for  turbulent  flows.  These  formulas 
are  only  available  if  ue  is  being  specified  (IOPT  = 1) . The 
basic  formulas  calculate  Cf  and  6 in  the  transform  plane 
(incompressible,  two-dimensional) . The  formulas  are  as  follows: 


Laminar  Flow 


Turbulent  Flow 


(50) 


(51) 


(52) 


(53) 
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These  formulas  provide  initial  values  for  boundary  layers  on 
flat  plates.  They  are  chosen  by  inputting  LSHAPE  = 1.  For 
other  geometries,  the  value  of  the  x coordinate  is  transformed. 
Thus,  LSHAPE  = 2 chooses  the  values  for  a circular  cylinder, 
where 


where  xa  is  the  axial  coordinate  and  LSHAPE  = 3 chooses 
the  values  for  a cone,  where 

X = Txa(Ti)  (55) 

These  formulas  have  been  found  to  be  quite  adequate  for  calcu- 
lating flows  over  long  bodies.  Small  initial  errors  in  the 
calculated  boundary  layer  become  negligible  in  a few  boundary- 
layer  thicknesses. 


8 . 0 PROGRAM  OUTPUT 

8.1  STANDARD  OUTPUT 

Several  options  are  available  for  output  from  the  programs. 
The  parameter  NPRINT  on  item  number  2 in  Table  II  chooses  either 
all  the  available  output  (NPRINT  = 1)  or  only  that  essential 
for  monitoring  the  progress  of  an  iterative  calculation.  The 
parameter  DXP  on  item  15  controls  the  printing  of  boundary-layer 
velocity  profiles.  For  iterative  calculations,  the  output  from 
the  inviscid-flow  program  is  controlled  by  NPRINT  only  on  the 
first  iteration  (Nl  = 1).  For  subsequent  iterations,  only  the 
short  form  output  is  printed  as  if  NPRINT  = 0. 

The  complete  program  output  is  presented  in  the  following 
list.  Steps  1-10  are  always  printed  for  the  first  iteration. 
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Of  the  remaining  steps,  those  denoted  by  an  asterisk  (*)  are 
those  printed  when  NPRINT  = 0. 

1.  Three-line  title  or  description.  . 

2.  List  of  all  values  of  integers  or  first  data  card. 

3.  List  of  Body  Geometry  input. 

4.  List  of  other  input  values  for  inviscid  flow. 

5.  List  of  input  values  for  boundary-layer  calculation. 

6.  List  of  Body  Shape  data  for  boundary  layer. 

7.  List  of  other  boundary-layer  input  quantities. 

8.  List  of  Variables  used  in  inviscid-flow  calculation. 

This  includes  the  body  geometry  and  other  parameters. 

9.  Computed  geometric  parameters  in  normal  direction  for 

inviscid  flow. 

J - Normal  grid  index. 

AN  - Normal  coordinate. 

G - Stretching  function  derivative  (refs.l  and  8). 

GH  - Stretching  function  derivative  at  half  intervals. 

10.  Computed  geometric  parameters  in  tangential  direction. 

I - Tangential  grid  index. 

S - Arc  length  along  reference  surface 
X - Axial  coordinate. 

Y - Radial  coordinate. 

THET  - Angle  of  reference  coordinate  surface,  0. 

For  closed  bodies,  9 is  the  same  as  the 
body  angle,  0B.  For  open  bodies,  0 = 0fi 
on  the  forebody  and  0=0  on  the  afterbody. 
THETB  - Body  angle,  0B. 

AK  - Surface  curvature  on  closed  bodies.  For  open 
bodies  AK  is  the  surface  curvature  on  the 
forebody  and  AK  = -(d2rw/dx2)  on  the  afterbody. 
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10.  F 
(cone. ) 


Derivative  of  the  tangential  stretch 
function  (refs.  1 and  8) . 


*11.  inviscid  relaxation  iteration  history. 

IT  - Iteration  number. 

DPMAX  - Maximum  0 correction,  max 

ij 

ID , JD  - I and  J location  of  DPMAX. 

rMAX  - Maximum  residual,  maxlR^jl,  where  R^j 

is  the  right-hand  side  of  the  difference 
equation. 

IR,JR  - I and  J location  of  RMAX. 

ISUB , ISUP  - Indicates  if  maximum  residual  occurred 
' at  a subsonic  or  supersonic  point. 


.IT  . IT— l 
0 . . — 0 . • 

11 


12. 


RAVG 

RFl 

QF3 

NS 

SEC/CY 


Average  value  of  the  residual. 

Relaxation  factor  for  subsonic  points. 
Relaxation  factor  for  supersonic  points. 
Number  of  supersonic  points. 

Time  for  iteration  cycle. 


List  of  solution  of  perturbation  potential. 


*13.  Tabulated  values  of  surface  pressure  coefficient,  Cp, 

Mach  number,  and  axial  velocity  on  the  body  along  with 
a rough  plot  of  Cp  along  the  body.  This  plot  is 
distorted  in  the  axial  direction  because  it  is  for 
equal  spacing  in  the  computational  space.  The  asterisks 
show  the  level  of  sonic  Cp. 

*14.  Drag  coefficient  by  trapezoidal  and  Simpson  integration 
of  the  cp * s on  the  real  body.  The  displacement 
surface  is  removed  for  calculation  of  the  drag. 
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15. 

16. 

17. 

18. 


*19. 


(NOTE:  If  the  two  values  differ  greatly,  it  is  probably 

due  either  to  the  computational  mesh  for  the 
inviscid  flow  being  too  coarse  or  the  inviscid 
relaxation  being  not  sufficiently  converged.) 

Mach  number  chart  of  the  flow  field  in  the  computational 
plane.  Numbers  printed  are  the  Mach  number  multiplied 
by  100.  The  I values  are  from  top  to  bottom  and  the 
J values  are  from  left  to  right. 

Coordinates  x and  y of  the  sonic  line. 

Boundary-layer  reference  velocity,  u_  , and  unit 

eo 

Reynolds  number,  Re  /L. 

co 

List  of  boundary-layer  quantities  with  profiles  at 


intervals 

governed  by  DXP. 

AX 

- 

Axial  distance  from  the  nose,  x. 

UTAU 

- 

Friction  velocity,  uT. 

DELTA 

- 

Boundary-layer  thickness,  6. 

DELST 

- 

Displacement  thickness,  <5* . 

THETA 

- 

Momentum  thickness,  9. 

CF 

- 

Skin-friction  coefficient,  Cf. 

UB/UZ 

- 

Free-stream  velocity  ratio,  ue/ue  . 

CP 

- 

Surface  pressure  coefficient,  Cp. 

DELST+R 

- 

Augmented  body  radius,  6*  + rw. 

DX 

- 

Integration  step  size. 

HTR 

- 

Transformed  shape  factor,  <5^/0^. 

Quantities 

showing  status  of  inner  iteration. 

XMAX 

Location  of  maximum  change  in  boundary 
layer  6*  from  previous  iteration. 

DPMAX 

- 

Maximum  change  in  6*. 

DSMAX 

— 

Value  of  6*  at  XMAX. 
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*20.  Quantities  needed  for  restart.  These  are  the  values 

required  on  input  card  number  4.  Note,  however,  that 
the  restart  procedure  requires  the  particular  set  of 
these  values  which  is  the  last  to  be  printed  before 
termination  if  the  restart  file  (data  file  5)  is  to 
be  used  (see  section  10.2). 

21.  Rough  plot  and  list  of  ue^  and  uei*  This  is  printed 
at  the  end  of  the  inner  iteration  cycle  whenever  a 
new  comparison  is  being  made  between  the  viscous  and 
inviscid  velocities. 

8.2  SPECIAL  OUTPUT  MESSAGES 

Several  special  messages  are  contained  in.  the  output  to  call 
attention  to  specific  conditions  that  may  occur.  The  messages 
are  listed  in  this  section  with  instructions  about  what  to  do 
when  they  are  encountered. 

(1)  — DIVERGENCE  . RMAX  EXCEEDS  RCHEK. 

GO  DIRECTLY  TO  JAIL.  DO  NOT  PASS  GO. 

DO  NOT  COLLECT  $200. 

This  message  is  printed  by  the  inviscid-flow  program  if  the 
relaxation  procedure  diverges.  Check  all  input  to  verify  that 
it  is  correct.  If  no  obvious  errors  appear,  the  difficulty  is 
probably  either  in  the  choice  of  parameters  for  the  computational 
mesh,  or  the  smoothness  of  the  data  defining  the  body  shape. 

(2)  RFl  DECREASED  TO  BECAUSE  10-CYCLE  AVG  FOR  RMAX  INCREASED. 

This  message  refers  to  the  subsonic  relaxation  factor  in 
the  inviscid-flow  program.  The  initial  value  is  1.4.  The  value 
is  automatically  reduced  by  10  percent  if:  (1)  the  maximum 

residual,  averaged  over  10  cycles,  is  greater  than  that  for  the 
previous  10  cycles  and  (2)  the  last  maximum  residual  occurred  at 
a subsonic  point. 
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(3)  QF3  INCREASED  TO  BECAUSE  10 -CYCLE  AVERAGE  OF  RMAX 

INCREASED. 

This  message  refers  to  the  supersonic  damping  factor  in  the 
inviscid-flow  program.  The  initial  value  is  0.1.  The  value  is 
automatically  increased  if:  (1)  the  maximum  correction,  averaged 

over  10  cycles,  is  greater  than  that  for  the  previous  10  cycles, 
and  (2)  the  last  maximum  residual  is  at  a supersonic  point. 

(4)  INPUT  FROM  TAPE13  HAS  INCOMPATIBLE  DIMENSIONS 

This  message  is  printed  if  the  dimensions  of  the' 
solution  read  from  Logical  unit  13  (data  file  2 in  figure  3) 
are  not  the  same  as  the  values  of  IMAX  and  JMAX  put  on  item 
number  6 in  Table  II. 

(5)  ****ITERATI0N  FOR  BOUNDARY  LAYER/INVISCID  FLOW  EQUILIBRIUM 

CONVERGED 

This  message  is  printed  whenever  the  maximum  change  in  6* 
between  iterations  is  less  than  the  specified  percent. 

(6)  METHOD  FOR  CALCULATING  UTAU  IN  DERIV  DOES  NOT  CONVERGE 

This  message  refers  to  the  iteration  used  to  solve  equa- 
tion (43)  for  UT  when  6*  is  prescribed  in  the  bound ary -layer 
calculat '.on.  The  only  known  cause  of  the  iteration  failing  to 
converge  is  an  error  in  the  input  data. 

(7)  DELTA I HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  PROFILE  AT  PREVIOUS  STEP 

This  message  refers  to  the  transformed  boundary-layer 
thickness,  6^.  The  error  condition  may  occur  due  to  the  initial 
integration  step  size  DXZ  being  too  large.  Another  possible 
cause  might  be  a too  sudden  change  in  the  body  shape,  or  in  the 
prescribed  ue  or  6*  distribution. 
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(8)  DELTA  HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  VALUES  AT  PREVIOUS  STEP 

This  message  is  not  expected  to  occur  in  the  finished 
program.  If  it  does,  check  the  input  data  carefully. 

(9)  INTEGRATION  INTERVAL  HAS  BECOME  TOO  SMALL. 

This  message  usually  occurs  when  the  singularity  at  UT  = 0 
is  approached  in  a positive  pressure  gradient  with  ue  prescribed. 
It  may  also  occur  under  certain  other  conditions. 

(10)  METHOD  FOR  CALCULATING  INITIAL  VALUE  OF  DELI  DOES  NOT 
CONVERGE 

When  initial  values  of  Cf,  6,  or  6*  are  known,  the  calcu- 
lation must  solve  an  integral  equation  for  the  initial  value  of 
the  transformed  thickness,  6j_.  This  is  done  by  iteration  in  a 
similar  manner  as  for  uT  described  in  message  (6) . If  the 
iteration  does  not.  converge,  it  is  usually  due  to  errors  in  the 
input  quantities. 

(11)  INTERMEDIATE  RESULTS  OF  THETA  ITERATION 

This  message  is  printed  at  the  end  of  each  step  of  the 
0S  cycle.  It  is  followed  by  the  current  value  of  the  0S 
increment,  DTHET  and  the  value  of  the  squared  deviation,  DIFFS. 

(12)  FINAL  RESULTS  OF  THETA  ITERATION 

This  message  is  printed  whenever  the  0S  cycle  discussed 
in  section  4.3  has  converged.  The  message  is  followed  by  a list 
of  the  values  of  XSEP,  DTHET,  and  the  squared  error  for  the  best 
solution  of  that  cycle. 

(13)  FINAL  RESULTS  OF  XSEP  ITERATION 

This  message  is  printed  whenever  the  least  squared  error 
has  been  found  for  the  xs  cycle  discussed  in  section  4.3. 

When  this  condition  occurs,  the  values  of  0,,  and  x are  set 
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to  those  corresponding  to  the  best  solution  and  the  calculation 
reenters  the  inner  cycle  to  perform  a final  iteration  sequence. 
The  final  output  of  the  inviscid— flow  and  boundary-layer  solu- 
tions then  correspond  to  the  final  converged  solution. 

(14)  **********SOLUTION  CONVERGED  TO  LEAST  SQUARED  ERROR  TOLERANCE 

This  message  is  printed  whenever  the  squared  error  is  less 
than  0.0005  at  any  stage  of  the  iteration.  The  inviscid  and 
viscous  solutions  immediately  preceding  this  message  are  then 
the  best  solutions  of  the  calculation  procedure. 

(15)  SKIN  FRICTION  HAS  BECOME  NEGATIVE  IN  AN  INCORRECT  MANNER. 

CHECK  ALL  INPUT  CAREFULLY 

This  message  will  be  printed  if  the  skin-friction  coeffi- 
cient changes  sign  suddenly.  This  has  been  known  to  occur  on 
rare  occasions  when  the  friction  was  approaching  zero  and  the 
singularity  did  not  cause  the  integration  step  size  to  decrease 
sufficiently  to  stop  the  calculation  before,  a negative  value  of 
Cf  was  calculated.  In  such  a case  the  calculation  may  be  all 
right.  The  condition  has  also  occurred  at  other  times  when  an 
input  error  caused  a discontinuity  in  the  prescribed  ue. 

9.0  PROGRAM  OPERATING  PROCEDURE 

In  this  section,  the  construction  of  card  decks  for  opera- 
tion of  the  computer  programs  is  described.  First,  a general 
description  of  the  operations  required  is  given.  Then  the 
specific  Job  Control  cards  needed  for  operation  on  an  IBM  370 
computer  are  listed.  The  same  card  decks  should  be  applicable 
at  any  370  installation  with  minor  modifications. 
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9 . 1 GENERAL  JOB  CONTROL  SEQUENCE 

The  following  list  is  the  general  Job  Control  procedure 
that  would  be  required  to  run  the  programs  for  a complete  viscid- 
inviscid  interaction  calculation.  The  reader  is  referred  to 
figure  3 and  Table  I. 

1.  Create  partitioned  data  sets  for  restart  files 

(files  1-5  in  figure  3) . 

2.  Define  units  2,  3,  8,  9,  and  10.  These  unit  numbers 

are  needed  for  output. 

3.  Define  units  11,  12,  13,  14,  and  15  if  NRSTRT  = 1 in 

the  input  data.  These  unit  numbers  correspond  to 
the  input  files.  They  contain  data  created  in  a 
previous  run. 

For  starting  an  initial  calculation,  the  partitioned  data 
sets  would  be  created  in  a separate  operation.  Then,  since  no 
data  would  be  on  file,  only  units  2,  3,  8,  9,  and  10  heed  to  be 
defined.  For  restarting  an  iterative  calculation,  all  data 
files  would  exist,  so  units  11  to  15  must  also  be  defined. 

To  execute  the  boundary-layer  program  alone,  unit  10  must 
be  defined  in  order  to  output  the  6*  + rw  and  ue/ueQ  list. 

Unit  12  must  be  defined  when  LVARl  = 2,  and  unit  15  must  be 
defined  when  LDSTAR  = 1. 

To  execute  the  inviscid  program  alone,  units  2 and  8 must 
be  defined.  Unit  13  is  also  required  when  LREADP  = 1,  and 
unit  14  is  required  when  IXY  = 0. 

9.2  JOB  CONTROL  EXAMPLES 

In  this  section,  specific  examples  of  Job  Control  cards 
used  for  the  operations  discussed  previously  are  presented.  In 
the  examples,  the  computer  program  is  referred  to  as  "ITER" 
with  the  source  code  names  "SITER"  and  the  load  module  or  binary 
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version  named  "BITER".  The  account  ID  used  in  the  examples  is 
WYL.XM.K01.  Logical  units  5 and  6 are  the  standard  input/ output 
file  numbers.  It  is  not  necessary  to  specifically  define  these 
unit  numbers  in  the  JCL  deck. 

9.2.1  Creating  Partitioned  Data  Sets 

partitioned  data  sets  for  use  as  input/output  disk  files 
must  be  created  before  the  normal  program  operation  can  proceed. 
The  following  procedure  is  suggested: 

Use  IBM  Utility  Program  IEFBR14. 

Use  default  values  for  DCB  (DSORG=PO, RECFM^VS) . 

On  3330  disk,  use  SPACE  in  tracks  as  follows  (refer  to 
figure  3 and  Table  I for  explanation  of  file  numbers) : 

SPACE  = (TRK, (2,1,10) ) 

SPACE  = (TRK, (10, 2, 10) ) 

SPACE  = (TRK, (20,4,10) ) 

SPACE  = (TRK, (4,1,10) ) 

SPACE  = (TRK, (6,1,10) ) 

Example  of  creating  a partitioned  data  set  called  VELBOD: 

//EXEC  PGM=IEFBR14 

//A  DD  DSN=WYL.XM. KOI. VELBOD, VOL*=volume, 

//  UNIT=3330,DISP= (,CATLG) , 

//  SPACED (TRK, (2,1,10))  , • . 

9.2.2  Starting  an  Iteration  Sequence 

To  start  an  iteration  sequence,  unit  numbers  2,  3,  8,  9, 
and  10  must  be  defined.  The  specific  sequence  of  cards  used  to 
perform  the  calculations  presented  in  section  10.1  is  as 
follows: 


VELBOD  (File  3) 
RESTRT  (File  5) 
PHI  (File  2) 
XOFILE  (File  1) 
DSFILE  (File  4) 
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//  EXEC  FORTGO, PROG= ITER, VOL^ volume 
//  LIB— 'WYL.XM.K01 .BITER  1 

//GO.FT02F001  DD  DSN=WYL.XM.K01 .VELBOD (RUNl) , 
//  DISP=OLD 

//GO.FT03F001  DD  DSN=  WYL.XM. KOI  .RES TRT  (RUNl) , 
//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL. XM. KOI. PHI (RUNl) , 

//  DISP=OLD 

//GO.FT09F001  DD  DSN=WYL .XM. KOI .XOFILE (RUNl) , 
//  DISP=OLD 

//GO. FT10F001  DD  DSN=WYL. XM. KOI .DSFILE (RUNl) , 
//  DISP=OLD 

//GO. SYS IN  DD  * 


Input  data  cards 


/■* 


9.2.3  Restarting  an  Iteration  Sequence 

The  specific  cards  used  to  perforin  a restart  of  the 
calculation  started  in  the  previous  section  are: 

//  EXEC  FORTGO, PROG= ITER, VOL= volume, 

//  LIB= 'WYL.XM. KOI. BITER* 

//GO.FT02F001  DD  DSN=WYL .XM. KOI .VELBOD (RUN2) , 

//  DISP=OLD 

//GO.FT03F001  DD  DSN=WYL.XM.K01.RESTRT(RUN2) , 

//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL .XM.K01 . PHI (RUN2) , 

//  DISP=OLD 

//GO.FT09F001  DD  DSN=WYL .XM. KOI .XOFILE (RUN2) , 

//  DISP=OLD 

//GO. FT10F001  DD  DSN=WYL .XM. KOI .DSFILE (RUN2) , 

//  DISP^OLD 
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//GO.FTllFOOl  DD  DSN=WYL  .XM.KOl  . RESTRT  (RUNl)  , 

//  DISP=OLD,LABEL=(,  , , IN) 

//GO.FT12F001  DD  DSN=WYL  .XM.KOl  .VELBOD  (RUNl)  , 

//  DISP=OLD,LABEL= {, , , IN) 

//GO.FT13F001  DD  DSN=WYL. XM.KOl . PHI (RUNl) , 

//  DISP=OLD,LABEL= ( , , , IN) 

//GO.FT14F001  DD  DSN=WYL .XM.KOl .XOFILE (RUNl) , 

//  DISP=OLD,LABEL=(,  , , IN) 

//GO.FT15F001  DD  DSN=WYL .XM.KOl .DSFILE (RUNl) , 

//  DISP=OLD,LABEL= (, , , IN) 

//GO. SYS IN  DD  * 

^ Input  data  cards 

/* 

These  cards  were  used  with  the  example  discussed  in  section  10.2. 

9-2.4  Executing  the  Boundary-Layer  Program  Alone 

The  specific  cards  used  to  perform  the  calculations  dis- 
cussed in  section  10.3  are  listed  in  this  section,  in  the 
example  shown  here,  all  input  is  assumed  to  be  from  cards,  but 

the  output  list  of  6* + r and  u /u  is  to  be  saved  on 
• w e “o 

unit  lo.  Unit  12  would  be  required  for  input  if  LVARl  = 2, 

and  unit  15  would  be  required  if  LDSTAR  = 1.  The  cards  used 

in  the  example  in  section  10.3  are: 

//  EXEC  FORTGO,PROG= ITER, VOL= volume 
//  LIB= 'WYL. XM.KOl .BITER' 

//GO.FTlOFOOl  DD  DSN=WYL .XM.KOl .DSFILE (RUNl) , 

//  DISP=OLD 

/ / GO . SYS IN  DD  * 


Input  data  cards 
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9.2.5  Executing  the  Inviscid  Program  Alone 


The  cards  used  to  perform  the  calculations  discussed  in 
section  10.4  are  listed  in  this  section.  In  this  example,  the 
velocity  potential,  0,  is  input  from  unit  13,  and  the  new  solu- 
tion for  0 is  output  on  unit  8.  The  calculated  velocity  on 
the  body  is  output  on  unit  2.  Input  from  unit  13  corresponds 
to  LREADP  =1  in  the  card  input  data. • In  addition,  unit  14 
would  be  required  for  input  of  the  body  shape  if  IXY  = 0 in 
the  card  input  data.  The  specific  cards  used  in  the  example  are: 


//  EXEC  FORTGO, PROG=  ITER, VOL= volume, 

//  L IB=  ' WYL  . XM . K0 1 . B ITI}R 1 

//GQ.FT02F001  DD  DSN=WYL.XM.K01 .VELBOD (RUN2) , 
//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL.XM.K01 . PHI  (RUN2)  , 

//  DISP=OLD 

//GO.FT13F001  DD  DSN=WYL.XM.K01 . PHI (RUNl) , 

//  DISP=OLD , LABEIr=  { , , , IN) 

//GO. SYS IN  DD  * 


Input  data  cards 


/' 


\ 

10 . 0 NUMERICAL  EXAMPLES 


In  this  section,  several  example  calculations  are  presented 
to  aid  in  program  checkout.  An  example  is  presented  of  a complete 
viscid-inviscid  interaction,  starting  from  the  initial  input  and 
restarting  after  premature  termination.  An  example  is  also  pre- 
sented of  the  use  of  the  boundary-layer  program  alone  for  a two- 
dimensional  geometry.  That  example  also  demonstrates  the  two 
options  for  boundary  conditions,  having  ue  specified  in  the 
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beginning  of  the  calculation,  and  6*  specified  in  the  second 
part.  Input  data  for  another  sample  case  are  also  presented  to 
demonstrate  the  use  of  the  program  to  calculate  the  inviscid 
flow  alone. 

10.1  AXISYMMETRIC  INTERACTION 

A list  of  the  punched  card  input  data  for  a sample  calcu- 
lation on  the  boattailed  body  shown  in  figure  4 is  presented  in 
figure  5.  The  case  being  calculated  is  for  a free-stream  Mach 
number  of  0.9.  The  body  corresponds  to  the  ogive-cylinder  body 
with  a circular-arc-conic  boattail  described  in  reference  11. 

The  JCL  card  deck  for  this  case  has  been  presented  in  section 
9.2.2.  The  running  time  for  the  complete  calculation  is  about 
10  minutes  on  the  IBM  370/165. 

Selected  output  for  the  sample  case  is  shown  in  figure  6. 

A total  of  27  iterations  were  calculated.  The  best  solution 
was  determined  to  be  that  corresponding  to  the  values  of  xg 
and  0g  used  for  iteration  20.  Those  values  are  used  in  itera- 
tion 27  to  produce  a final  list  of  the  solution  so  that  it  is 
not  necessary  to  search  through  the  output  for  the  best  solution. 
However,  some  caution  is  required  in  interpreting  the  final  out- 
put. The  agreement  between  the  ue^  and  the  ue^  for  the  final 
inviscid  and  viscous  solutions  is  somewhat  dependent  upon  the 
accuracy  of  the  previous  calculations  as  determined  by  the  input 
quantities  NlMAX  and  DFACT . For  the  calculations  shown  in  fig- 
ure 6,  NlMAX  had  a value  of  7 and  DFACT  had  a value  of  0.5.  A 
smaller  value  of  NlMAX  and  a larger  value  of  DFACT  should  improve 
the  final  agreement. 

The  complete  list  of  output  for  this  case,  with  NPRINT  = 0, 
consisted  of  a total  of  6661  lines.  Output  steps  1-7  as  listed 
in  section  8.1  have  been  omitted  from  this  presentation  since 
they  simply  verify  the  input  data.  The  output  pages  shown  are 
those  corresponding  to  steps  8-21  of  the  set  described  in 
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section  8. If or  NPRINT  =0  on  the  first  iteration  and  then  a 
few  samples  of  the  rough  plots  of  and  uei  at  intermediate 

steps,  including  that  at  iteration  20  (fig.  6 ( j ) ) . Finally, 
figure  6 concludes  with  the  inviscid  solution  and  boundary-layer 
solution  corresponding  to  the  final  result. 

10.2  EXAMPLE  OF  RESTARTING  AN 

INTERACTION  CALCULATION 

The  punched  card  input  data  for  restarting  the  interaction 
calculation  of  section  10.1  after  premature  termination  at  the 
twelfth  iteration  (eleven  complete  iterations)  is  presented  in 
figure 07.  . See  figure  6(i)  for  the  printed  output  of  these  quanti 
ties.  The  JCL  card  deck  for  this  calculation  was  presented  in 
section  9.2.3.  Note  that  the  iterative  calculations  cannot  be 
restarted  at  any  arbitrary  iteration  using  the  restart  file, 

unit  11.  That  file  and  the  other  output  files  contain  only  the 

0*  *.  * 

data  that  were  output  just  prior  to  the  termination. 

10.3  TWO-DIMENSIONAL  BOUNDARY  LAYER. 

A list  of  the  punched  card  input  data  for  a sample  calcu- 
lation 'on  the  two-dimensional  configuration  shown  in  figure  8 
is  presented  in  figure  9.  Note  that  two  sets  of  input  are 
presented  for  this  case,  giving  an  example  of  the  options  of 
prescribed  ue  and  prescribed  6*.  The  data  for  u0  and  6* 
were  obtained  from  the  experimental  results  of  reference  12 
which  indicate  separation  occurring  in  an  adverse  pressure 
gradient  region  downstream  of  a shock  wave.  The  output  for  the 
complete  boundary-layer  calculation  are  presented  in  figure  10. 

No  external  data  files  were  used  for  input  for  this  case.  The 
JCL  card  deck  for  this  case  was  presented  in  section  9.2.4. 
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10.4  AXISYMMETRIC  INVISCID  FLOW 

The  punched  card  input  data  for  a sample  calculation  of 
the  inviscid  flow  alone  are  presented  in  figure  11.  These  data 
correspond  to  the  first  inviscid-flow  step  of  the  iteration 
calculation  described  in  section  10.1.  The  output  for  this 
case  is  the  same  as  that  shown  in  figure  6(a)  through  (e) . 
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TABLE  I.  RELATION  BETWEEN  EXTERNAL  DATA  SETS 
AND  INPUT/OUTPUT  LOGICAL  UNIT  NUMBERS 


Data  File 
(fig.  3) 

Logical  Unit 

Input 

Output 

1 

14 

9 

2 

13 

8 

3 

12 

2 

4 

15 

10 

5 

11 

3 
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TABLE  II.  INPUT  DATA  CARDS 


CARD 

ORDER 

VARIABLES 

FORMAT 

1 

TITLE  (three  cards) 

20A4 

2 

LPROG , NRSTRT,  NPRINT , NlMAX , Nl.N2.N3, IBL , MIT2 

1615 

3 . 

GAMMA, AMINF 

8F10.0 

If  LPROG=l,  skip  to  item  no.  11. 

If  LPROG=0,  put  in  all  cards  as  required. 

If  LPROG=0  and  NRSTRT>0,  stop.  No  further 
cards  are  needed. 

4 

IXY 

1615 

If  IXY=0,  skip  next  card. 

5 

X0,Y0 

2F10.0 

6 

IMAX , JMAX , MIT , MHALF , KLOSE , NPLOT , LREADP 

1615 

7 

DETAO.XIXM.XM,  XIMl  ,XBT  i 

8F10.0 

If  LPR0G=-1,  skip  remaining  cards.  , 

8 

XZNEW , DTHET 

8F10.0 

9 

LSEP 

L5 

If  LSEP= FALSE,  skip  next  card. 

10 

XSEP 

8F10.0 

11 

IOPT , K , LVARl , LSHAPE , L IC , LDSTAR , IUNIT , LSHPBL 

1615 

If  LPROG^l  or  if  LPROG=l  and  both  I0PT=1 
LVARl=2,  skip  items  12  and  13. 

12 

NVAR 

1615 
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TABLE  II.  CONCLUDED 


CARD 

ORDER 

VARIABLES 

FORMAT 

13;' 

XVAR,  VAR 

8F10.0 

14  . 

EL, PT, TT, TWONTT, VISC, RGAS , SCON, DFACT 

8F10.0 

15 

XZ , RLEN , XT , DXZ , DXP , DXOUT , HL IM 

8F10.0 

If  LSHAPEj^O,  skip  next  card. 

If  L IC= 1 , input  CFCl  and  DELTAl. 
If  LIC=2,  input  CFCl  and  DELSTl. 

16 

CFCl, DELTAl  (or  DELSTl) 

8F10.0 

If  IOPT^l,  skip  next  card. 

■ 

17 

UEl,DUEDX 

8F10.0 

If  LSHPBL=0  and  LPROG=0,  skip  items  18  and  19. 

* 

18 

NR  ' i 

1615 

19 

XRP,  RL 

■ 

2F10.0 
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Figure  1.  Boundary-layer  coordinate  system. 


64 


Figure  2.  Schematic  of  iteration  procedure. 


65 


AEDC-TR-77-72 


Figure  3.  General  relationship  of  programs  and  data  files. 
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' 'T* 7.68  > 


(b)  Detailed  boattail  geometry. 


Figure  4.  Axisymmetric  body  for  sample  calculation. 


67 


AEDC  TR-77-72 


BLAHA, CHAMBERLIN  AND  BOBER  6524  BOATTAIL-STING 
TESTING  COMPLETE  ITERATION 

TRANSONIC  FLOW  WITH  WEAK  SHOCK  ON  AFTERBODY 

0 0 0 7 1 1 1 1 „ 10 


1.4 

0.9 

53 

0.0 

0.0 

0.25- 

6.085 

0.50 

0.17 

0.75 

0.255 

1.0 

0.34 

1.5 

0.49 

2.0 

0.64 

4.0 

1.2 

6.0 

1.68 

8.0 

2.13 

10.0 

2.56 

12.0 

2.92 

14.0 

3.24 

16.0 

3.5 

18.0 

3.65 

20.0 

3.80 

22.0 

3.91 

24.0 

4.0 

26.0 

4.0 

30.0 

4.0 

35.0 

4.0 

40.0 

4.0 

45.0 

4.0 

50.0 

4.0 

50.4 

4.0 

50.5 

3.9997 

50.7 

3.997 

51.0 

3.988 

51.25 

3.9759 

51.5 

3.9597 

51.75 

3.9392 

52.0  . 

3.9146 

52.5 

3.8525 

53.0 

3.7733 

53.5 

3.6767 

54.0 

3.5624 

54.5 

3.4298 

(a)  First  43  cards. 

Figure  5.  Input  data  for  interaction  calculation 

on  body  of  figure  4. 
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55.0 

3.2786 

55.5 

3.1080 

56.0 

2.9175 

56.5 

2.7061 

57.0 

2.4835 

57.5 

2.2609 

58.08 

2.0 

58.5 

2.0 

59.0 

2.0 

59.5 

2.0 

60.0 

2.0 

60.5 

2.0 

61.0 

2.0 

62.0 

2.0 

63.0 

2.0 

64.0 

2.0 

51  31 

10  0 

0 

0.25 

0.84. 

58.08 

52.0 

0.0 

F 

1 1 

2 3 

1 

1.0 

16.81 

590.0 

2.0 

64.0 

6.0 

0 

80.0  58.08 


2 ; 0 
0.97 

0.001  0.0  0.5 


(b)  Remaining  23  cards. 
Figure  5.  Concluded. 
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VALUES  OF  VARIABLES  USED  IN  I NV I SCI D ‘ FLO*  CALCULATION 


O 


X FAX  JM AX  WIT  MHALF  KLObL  NPLOT  R:A OP 


SI  31 

20  U 

COO 

RF13 

0 • 1 AO  £♦01 

COVERG 
O.lOOc  + Ol 

OF  3 

0.100c*00 

OSDXIO 
C • 1 OOE+0  2 

DNDZG 
0 «S74E+0 I 

XXXM  XM  3SDX1M 

0 • 64QE+00  0.581E*02  0. 1915*02 

GAM 

0 • I AO  E + 01 

AMINF 

0.900E+00 

kstar  PSTAR 

CPSTAR 

CP  J 

46  0«893SE*00  -0 • 1 ti 79 c + 00 

0.l219ct0l 

RMSQ 

AMSQ 

GM1 

0 • 1 60  OE*  02  0 • 81 00E+0D 

0*40  00&+00 

oxio 

DZ 

DX  I SO 

0 • 2000c- 0 1 0 • 3333c—  01  0.2500E+04 


0M1U2 

0 #200  Qc  *00 
DX  10  2 

0*  7S00E  + 03 


GOGH l GMSQ 

0.3500E+01  0*1 620c ^ 00 

OZSQ  DZ2 

0 • 90  00E+0  3 0*1 500E+  02 


AOSO 

0«14352*0I 
DX  12 

0.2500E+0 2 


(a)  List  of  variables  used  in  inviscid-f low  calculation. 

Figure  6.  Selected  output  for  calculation  of  viscid- inviscid 
interaction  on  boattailed  body. 
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N0PM4L  COOBD.  STOETCH  FOB  ALFs  1.300 


J 

AN 

G 

Gh 

1 

—0.1 36  IE 

22 

o.? 

r .27*42— C4 

' ""  * 

4 

r 

03 

'i.54  «*E-C4 

0.1  ft 2E-03 

3 

c . ; 5 1 1 r 

03 

C .26845-C3 

C .4  7795  — C 3 

4 

C.1Q31E 

03 

C .6*745-03 

C.l 0155-C2 

5 

C.683CE 

32 

C.l  143S-C2 

0.1S-2E-C2 

' 

6 

C .49^45 

r2 

.22615-' 2 

C .25645-02 

7 

C • 3722E 

.*» > 

0 . 34675*02 

f .42 255-C2 

. - - 

a 

C .2919? 

0 2 

C .4 9&25-C2 

T.5SC6E-G2 

9 

C .3347E 

0 2 

C .fc8295-<’2 

C.7928E-C2 

10 

_0  • ? 9.2  2«t 

y? 

0. 90295-02 

o.l 3315-01 

1 1 

C • 1 597E 

C 2 

0.1 1695-01 

0.1  2:05-0 

12 

_C  • 1 ’40E. 

02 

A.1456E-fl 

0. 1 62?  5-**l 

13 

0.1!  34E 

02 

C.  179 45- Cl 

0.1 98 5 E— 01 

1 4 

C . ^6  4 9E 

01 

o.  21 75E-C1 

0. 23885-C1 

% 

15 

C . 4J47F 

Cl 

C .26'’  25-01 

0.2  6395—0 1 

16 

_ C •▼.:6 9c 

01 

V.39  76*:-C1 

0 • 33375— Cl 

17 

C .60666: 

01 

0 . 259  95-01 

0.  J847E-01 

ie 

0 • 5.?  . 65 

:i 

C .41 745-01 

0 .44895— ri 

19 

T.4aa!F. 

ci 

0. 48*35-01 

C .51*65-01 

20 

C . 3b  1 2E 

Cl 

C • 54  r 95  — C 1 

0 .58605-C1 

21 

0 • 12425 

:i 

0. 62215-0 1 

0.66325-0 1 

2 2 

C .27346. 

31 

r .70  355-0 1 

C.74AQE-CI 

23 

C.2291E 

bt 

0.79C  25  — 0 1 

0. 8369E-C1 

24 

__i  • 1.8926 

? .283  55-0 1 

0.93355-0! 

25 

C .1535c 

■51 

C .98 15E-C1 

C.1C375  CC 

2o 
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(d)  Iteration  history  for  inviscid  relaxation  on  plain  body. 
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Figure  6.  Continued. 
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(g)  End  of  first  boundary- layer  calculation. 
Figure  6.  Continued.  . 
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Figure  6.  Continued. 
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(i)  Comparison  of  viscous  and  inviscid  velocities  for  an 
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Figure  6.  Continued 
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Figure  6.  Continued. 
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Figure  6.  Continued. 
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Figure  6.  Concluded, 
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Figure  7.-  Input  data  for  restarting  calculation  of  viscid-inviscid 

interaction  after  4 iterations. 
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15480. 

10.0 

16100. 

10.25 

14800. 

10.75 

13400. 

11.0 

12600. 

11.50 

12100. 

12.0 

11150. 

13.0 

10900. 

14.0 

10800. 

15.0 

10300. 

16.0 

9900. 

17.0 

9800. 

18.0 

10200. 

1.0 

15.0 

0.0 

15.0 

.00233 

.065 

19 

0.0 

0.0 

1.0 

0.0 

2.0 

0.0 

3.0 

0.0 

4.0 

0.25 

5.0 

0.50 

6.0 

0.70 

7.0 

0.85 

8.0 

0.98 

9.0 

1.0 

10.0 

0.98 

11.0 

0.85 

12.0 

0.70 

13.0 

0.50 

14.0 

0.25 

15.0 

0.0 

16.0 

0.0 

17.0 

0.0 

18.0 

0.0 

2 0 


585.0 

-5.0 


2 0 


.99  0.0 

001  0.0 


0.0 

0.5 


0.0 


0.0 


(a)  Input  for  case  with  ue  specified. 

Figure  9.  Input  data  for  two-dimensional 
boundary-layer  calculation. 
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ALBER  TEST  CASE 

TWO-DIMENSIONAL  BOUNDARY  LAYER  CALCULATION 
SPECIFYING  DELST  DISTRIBUTION  FROM  X=10.75 

1 

1.4  .7325 

20002020 

13 

10.0  .03 

10.25  .05 

10.75  .07 

11.0  .083 

11.5  .112 

12.0  .148 

13.0  .314 

14.0  .435 

14.5  .460 

15.0  .440 

16.0  .306 

17.0  .230 

18.0  .190 

1.0  15.0  585.0  .99  0.0  0.0  0.0  0.0 

10.75  18.0  -5.0  .001  0.0  0.5 

.00077 

13400.  -2820. 

19 

0.0  0.0 

1.0  0.0 

2.0  0.0 

3.0  0.0 

4.0  0.25 

5.0  0.50 

6.0  0.70 

7.0  0.85 

8.0  0.98 

9.0  1.0 

10.0  0.98 

11.0  0.85 

12.0  0.70 

13.0  0.50 

14.0  0.25 

15.0  0.0 

16.0  0.0 

17.0  0.0 

18.0  0.0 


(b)  input  data  for  case  with  6*  specified. 
Figure  9.  Concluded. 
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UZ  * 9.90272E  C3 

REZ  s 3.0  2426E  C5 


AX 


UTAU 


oTC  3T1X  24  2 = 0 2 

J5.C7999E-01 2*  21 451  E 02 

1 iCXCOOE  THJ  3 • 1 0 8 97E  ~0  2 
1 «5959Sg  00  3.0eci4E  C2 

2 »C  9 "»?  4E~00  2 • 990  19E '"OX” 

2.627^3=  00  2.59925E  02 

3.14  789=  00  2.f3  22eE“C  2“ 

3>yt8g  00  2.45522F  02 

4.3?764E~00  2.807905  C2 

_4.pq783E  00  3.305431  C2 

5.a:37  9E_  00~  3. 943756  02 
5.97^7 gp  00  4.64449E  02 


00 

<X> 


6.4C»77£  00 

_7.00  376E  00 
7.51 5765  00“ 
_fl  .C 3 5 7 2 E 0 0 
fl .5477  2E  00“ 
_9 .0  5971  E 00 
9.3717  5f  00 
1 .007575  oi 


5.07260"  C2 
5.46295c  02 
5.82  520  E 02 
6.19694E  C 2 
6.64  967-  "02“ 
7 .07950  E 02 


7.2e*3  <JE  02 
7.15222*  02 


l .062765 

01 

5.44731= 

C2 

2 

1 . 1 IA76E 

01 

4.23246= 

C 2 

2 

1.16596c 

01 

3.80  214E 

02 

■a 

1 .21 796E 

0 1 

2.59 195E 

02 

3 

1 .27555E 

0 1 

2.84449E 

02 

3 

1 » 3 2c7  55 

01 

2. 04469E 

02 

3 

1.28435= 

01 

a.adjesE 

02 

3 

! .43^156 

01 

2.72669E 

C 2 

4 

468355 

0 1 

2.51832E 

02 

4 

1 .501 15E 

01 

2.46391E 

02 

4 

DELTA 

3#  6 9 e tce-TT 

4.C02C4E-C1 

4.1  1 624E-C1 
4.22195E-CI 

4. 2 2d94E- Cl 
4.  t296£c-01 
4.33C4CE- Cl 
4.99273F-01 
"4,e«;e4  3F- ci 
4.7€7e?e- ci 
4. 76346E-CI 
4.9P0S6F-01 

“ 5 . 26  8 9 5E— ~G  t ~ 
5.73291E-01 
6.6 15  14E-01 
3.60215E-C1 
3.7764EE- Cl " 
2 .963  lcE- 01 

5i9a213C-oT“ 

3.57844E-01 


►15252E- Cl 


DEL ST 

6.500115-02 
6.SC991E-C2 
7. 34969E-C2 
7.555C1E-02 
7. 97424E-02 
_9. 921 06F-02 
1.14  8 C6E-6T* 
1 • 09046E-01 
9.  I 9850c -02 
7. 2C426E-02 
5.62S93E-02 
4. 4c252E-02 
4. 012C7E-02 
2. 72625E-02 
3.561 CCE-02 
7.44262C-02 
3. 3271 Oc  — 02 
3.21225E-0 2 
3. 26993E— 02 
3.421 83E-02 
4.6027JE-02 
6.24852E-C2 
7. 60te9E-0 2 
9. 62234E-C2 
1.C2443E-C1  " 
1 .04649F-GI 
1 .05517E-01 
1.13605E-CI 
1.26391E-01 
1.29893E-01 


THETA 

3.91  983  E— 02- 
4.  15509E-02 
4. 40558E— 02 
4.532  85  E —0  2 
4.77122E-C2 
5.70932E-02 
6.  AT 94  74  E -02“ 
6.30893E-02 
5.52  31 9E  — 02 
4. 493  03E-02 
3.43773E-02 
2.59309E-02 
2.21274E-02 
1 .92269E-02 
1 • 6 7643E  — 02 
1 • 72O59E-02 
1 .57975E-02 
_l .437166-02 
1 .43427E-02“ 
1 .54963E-02 
2.35717E-02 
3.27317E-02 
3. 90455E-02 
4.79382E-02 
5. 03771 E— 02 " 
5.20264E-02 
5. 32714E-02 
5. 69 322 E— 02 
6. 21 547E-02 
6. 35316 E “0  2 


CF 

2. 3'3d06E^O"3“ 
2.25014E-03 
2.171  I 2 £-0 3 

2.  153946-03 
2.03831 E— 0 3 
1 • 76648E-03 
1.540  9 2 6-0  3“ 
1.67291E-03 

*1 .9  9671E-0  3 
2.37  327E-0  3 
2.7644 3E-03 

3 .0  70  23  E-0  3 

3.1  75  17  E-0  3 

3.  274C5E-0  3 
3.3i934E-03 
3 . 349  29  E-0  3 
3.43576E-03 
3.5  12  56  £-0  3 
3.523776-0  3 
3.466416-03 
2.82507E-0  3 
2.  1o594£-03 
1.77953E-03 
1 .31879E-0  3 
1.31  428E-C  3 
1 .338716-03 
1.39131  E-0  3 
1 .30801  E-03 
1 • 1 70  1 1 E-0  3 
l .135446-03 


__  UE/UZ 
0 .96458 
•0  .97C76 
0.95715 

0 .95250 
* 0 .94026 

0.89244 
0.66059 
0.86980 
0 .906C8 
0.97186 

1 .0  6 362 
l .17270 
1 • 24  3<»5 
1 .31 123 
1 .37522 
1 .4J573 
1 • 50  590 
1 .56  70 1 
1*. 59929 
1 .58716 
1.38871 
1.25759 
1.19166 
1.12144 
1 . 10690 
1 .09801 
1 .09219 
1 .07196 
1 .04605 
1 .03966 


CP 

0.0  30  73 
0.05802 
O.Ce4  74 
0.09389 
0.11767 
0.20822 
0.26072 
0.20864 
0.14172 
0 .0 1524 
-0 . 1 5989 
-0.391 13 
-C.53C33 
-0.67677 
-0 • 8C889 
-0.94256 
- 1 .06694 
-1  • 19331 
-1 .25931 
-1.23561 
-0.83271 
-0.54985 
-0.41347 
-C. 26477 
-0.25584 
-0.24692 
-0.24879 
-0.19934 

-0.1 4oeo 

-0.11762 


OEL  ST  +R 

DX 

H TP 

6*500 FIE—  02 

i.0CC5-03~~ 

I.Tt 24“ 

6.90SS1E-02 

6.  400E-02 

1.2592 

“ 7. 34S69E-02 

6. 40C5-02 

l • 3 Co2 

_ 7.55SC1E-03 

1 .2  60E- 01 

1.3074 

7 • 974  24E—  02 

1 • 600E-02  ~ 

1.3218 

9. 821 C6E- 02 

l .28CE-C1 

1.4095 

iiEl7eCE-dl” 

3.  2 C CE  - 02 

1.4  796“ 

2.9801SE-01 

1 .280E-01 

1.4314 

4.18945E-CI 

6.4CCE-  32 

1.3355 

5.44COOE-01 

l .28CE-01 

1.2341 

t.3704fc£-01 

6. 4 COE  — 02 

1.1383 

7.4OStlE-0l 

1 .26CE-01 

1 • 0539 

a.  1 seetE-o: 

1 . 28CF-C1  ~ 

1.C071 

8. 877S3E- 01 

1.26CE-C1 

0.961 3 

S.5265SE-01 

1.28CE-C1 

0.9C5* 

1*0151 4E  00 

6.4C0E-C2 

1 .C5C7 

1.02423E  00 

6. 400E-C2 

1.0231 

1.C3CS3E  00 

1.280E-01 

0.9961 

1.02127E  00 

’ 6.40CL-02 

0.9990 

1.0043EE  00 

6. 0CCE-C3 

1.0344 

9.44<.34E-01 

6.4  0 Ofc—  C2 

1.1  95  9 

8.9034SE-01 

3.2C0E-C2 

1.3431 

8.271’IE-Cl 

1 . oCGE-02 

1.4450 

7.623136-01 

3.2C0E-02 

1 .6009 

e.E133eE-0t 

l .28CE-C1 

1 .5963 

5*  377C9E— 01 

O.40CE-C2 

1 .5803 

3.94639E— 01 

1 .28CE-C1 

1.5542 

2.70E17E-01 

6. 4 COE- 02 

1.5829 

1 .55525E-01 

1 .28CE-01 

1.6385 

1 . 29893E- 01 

6. 4 COE -02 

1.0537 

Figure  10. 


(a)  Output  from  case  with  ufi  specified. 

Output  from  two-dimensional  boundary— layer  calculation. 
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uz 

REZ 


= 9.9^2729  C 3 
±^3»y2*26E  v*5 


AX 

utau 

Cc-TA 

CELST 

TH£  T A 

CF 

UE/UZ  CP 

CELST+9 

HT9 

: e 

■?  1 

2.754485 

•*  2 

1 . ▼6*31 ?-Cl 

7.0. *'  f E-r2 

2 .579295-52 

7.6S997E-C4 

1 . 353! 6 -r .76251 

9.52  5*' 

!.(  *'E  - * 3 

2.1  ’ 3 

1 • : 259  £ 

”*  • 

• 

! • 7 Z4E  2 £ 

_ 2 

r.’ 

9 .e  . 2*  1E-C2 

2.1  75785-02 

3.41  *»S3E-l  4 

1.2853?  -C  -*.  *■»* 

9.C91 75F-02 

3.2C -£-r  ? 

2.4*38 

l •! 77\^£ 

v'l 

1 . .27385 

( 2 

2 • 7 ’ 1 iec-0  3 

1 • 2 5*  9t-«  1 

3.6802'.  E-02 

1 .3227Jfc-v4 

1.23465  -'.'.49*54 

“.6586*5-0! 

6 .4*-.*  =-  •? 

2.9*61 

1 • 

'■  » 

1 .‘  5375c 

• l 

3.0*M  35-01 

1 .*76225-61 

4..?f.7r45— 02 

1 .5*.??SE-^6 

1.19521  -r.  .4*7?4 

6. 37636E-C1 

3.2-  i E- * 2 

3.89^6 

UHl  9C 

’*1 

-6 . 1 £976  5 

C 1 

4 • 72E-C 1 

2.6ZM2E-'  1 

4 . 5Hw4~E  — C 2 

-5. 337*5F-'  5 

1 • 16732  -:'.37aji 

b.?C42?p-: i 

* .4?“ 9-' 2 

6. 3 ‘ *8 

1 . 7322  “=• 

o : 

< i 

5.9A-.  93E-01 

3.5?“ 71  £ -Cl 

4.729I2E-C2 

-1.1! 729E-C4 

1 . ! 74^6  -C  .369.'  7 

7.72346F-.'  1 

fc.4 rrc-^2 

t.54-1 

i .3  6-8  96 

j: 

-1  • . i 7c  riE 

0 2 

t.r-36  ♦?£-?! 

4.22766E-01 

4.C95345-C2 

-1 .491 61E-C4 

3.1667C  -C. 361  75 

6.9  3f'4?E-iVv 

l .26-9- “1 

7 # 52  6 2 

1.-4! '?£ 

«>: 

-4.il *5!  5 

vl 

7 * «*  1 3 .EE'-cl 

4.E5544F-01 

5. 4fc862E— C 2 

-1  .421  33E-f  4 

1.1544c  — i.  .3545’ 

6.*  26  2*5-''  • 

6 .4<  r 3 

7.4*  79 

1 .A3J29E 

01 

-S.lA2245 

• 1 

7. ““AfJE-Ci 

4 . A JC  erE-Cl 

6.091 4 25-02 

-l  ••'las.e-?* 

1 .12233  -C .’4623 

4 . c 2 .?  7 1 F — *.»  1 

e.4  ?rp-.'2 

6 • 4 * ’ 5 

1.  54  *4  9£ 

Ci 

-3 .8 ! 51 2 E 

.*! 

6.7E.  3FE-C1 

3 . 91 732t-?l 

6.  -973^5 -J2 

— 2 • 2?3  7 ) E— T 5 

1.1729!  -C • ? 35?  4 

3.81772F-?! 

6.4<V*s-'2 

4.7^74 

1 .534685 

: i 

3.5  5fc97  £ 

v 1 

< • ! 6 3925-01 

2.121 ?5E-r  1 

7.62d2dc-02 

2. 31262E-' 5 

! • “ 766  8 -0  . 1714a 

3 . 1 3! 255-C 1 ’ 

6.4-r  c.*2 

3.45-M 

1 • ->4«3.;c£ 

•m 

fc • 66549c 

r « 
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8.26252E-C 2 
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: .0f,?7*  -*  .? 

2.6671 7 « 

6.4?^F-'2 

2.7776 

l • TC 2SE 

01 

i • .3  7 v 1 4 ri 

0 2 

5.6?  - 5 8c-C ! 

2.Z99C2E-0; 

8.2535SE-02 

3 . 4^71 9E— ?4 

1 • 0549?  -C  • i 2697 

2.299525-01 

?.2cr.F--'? 

2.  93 

1*  y5';21  E 

Cl 

1 .7 1CT2£ 

2 

5 . 552555— •-  «. 

2 • “7517E-vl 

ti. 5761 6£— c2 

5. 369865-.*  4 

!.'S95c  -6.122?? 

’.n7517C-' | 

6 .4  -v  F--2 

2 • ? 52  * 

1 .3r  *%  :e 

c; 

1 75»F 

2 

5.5.'‘c5tE-C: 

1 .MftCCC-Cl 

8.23511E-92 
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* .94?A 

00 

-vl 


(b)  Output  from  case  with  specified. 
Figure  10.  Concluded. 
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BLAHA, CHAMBERLIN  AND  BOBER  6524  BOATTAIL -STING 
INVISCID  FLOW  ONLY 

TRANSONIC  FLOW  WITH  WEAK  SHOCK  ON  AFTERBODY 
-10  0 


1.4 

53 

0.9 

0.0 

0.0 

0.25 

0.085 

0.50 

0.17 

0.75 

0.255 

1.0 

0.34 

1.5 

0.49 

2.0 

0.64 

4.0 

1.2 

6.0 

1.68 

8.0 

2.13 

10.0 

2.56 

12.0 

2.92 

14.0 

3.24 

16.0 

3.5 

18.0 

3.65 

20.0 

3.80 

22.0 

3.91 

24.0 

4.0 

26.0 

4.0 

30.0 

4.0 

35.0 

4.0 

40.0 

4.0 

45.0 

4.0 

50.0 

4.0 

50.4 

4.0 

50.5 

3.9997 

50.7 

3.997 

51.0 

3.988 

51.25 

3.9759 

51.5 

3.9597 

51.75 

3.9392 

52.0 

3.9146 

52.5 

3.8525 

53.0 

3.7733 

53.5 

3.6767 

54.0 

3.5624 

54.5 

3.4298 

(a)  First  43  cards. 

Figure  11.  Input  data  for  inviscid-flow  calculation 

on  body  of  figure  4. 
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55.0 

3.2786 

55.5 

3.1080 

56.0 

2.9175 

56.5 

2.7061 

57.0 

2.4835 

57.5 

2.2609 

58.08 

2.0 

58.5 

2.0 

59.0 

2.0 

59.5 

2.0 

60.0 

2.0 

60.5 

2.0 

61.0 

2.0 

62.0 

2.0 

63.0 

2.0 

64.0 

2.0 

51  31 

10  0 

0.25 

0.84 

0 0 
58.08 


0 

80.0 


(b)  Remaining  18  cards 
Figure  11.  Concluded 
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NOMENCLATURE 

A factor  defined  by  equation  (13) 

A^.  coefficients  in  equations  (38) , (39)  and  (40)  to  (42) 

a speed  of  sound 

B factor  defined  by  equation  (14) 

C Chapman-Rubes in  parameter,  equations  (22)  and  (27) 

Cf  skin-friction  coefficient 

Cp  surface  pressure  coefficient 

H total  enthalpy 

Htr  transformed  bound ary- layer  shape  factor,  equation  (44) 

k factor  denoting  two-dimensional  (k  = 0)  or  axisymmetric 

(k  = 1)  flow 

L reference  length 

M Mach  number 

P Prandtl  number 

p pressure 

r radius 

S total  enthalpy  parameter,  equation  (4) 

s squared  error  factor,  equation  (49) 

T temperature 

Ts  factor  in  Sutherlands  viscosity  relation,  equation  (28) 

t axisymmetric  transformation  factor,  equation  (15) 

U, V transformed  velocity  components,  equations  (18)  and  (19) 
Up  wake  velocity  factor,  equation  (32) 
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U 

T 

U,V 


u,v 
X,  Y 
x,y 


x,y 


Y 


4- 


a 

P 

7 


* 


6 


★ 

k 


A 

P 

V 

p 


NOMENCLATURE  (CONTINUED) 
friction  velocity,  equation  (33) 

velocity  components  of  physical  flow  field,  figure  1 

transformed  velocity  components,  equations  (9)  and  (10) 
transformed  coordinates,  equation  (17) 
coordinates  of  physical  flow  field,  figure  1 
axial  location  of  peak  pressure  downstream  of  boattail 

axial  location  of  separation  point 


transformed  coordinates,  equations  (5)  and  (6) 


coordinate  of  logarithmic  part  of  boundary  layer, 
equation  (34) 

damping  factor  for  viscid-inviscid  iterations,  equation  (48) 
eddy -viscosity  factor,  equations  (35)  to  (37) 
ratio  of  specific  heats 
boundary-layer  thickness 


boundary-layer  displacement  thickness , 


1 . -£2L 

P " 


u \ r 

u J r 
e/  w 


dy 


boundary-layer  displacement  thickness 


■ /(■  - *) 


dy 


angle  of  body  surface  with  the  axis;  also  boundary-layer 
momentum  thickness 

constant  in  Sutherlands  viscosity  law,  equation  (28) 
molecular  viscosity 

W] p 

density 
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NOMENCLATURE  (CONCLUDED) 


t shear  stress 

stream  function 

SUBSCRIPTS 

e refers  to  bound ary- layer  edge 

I refers  to  inviscid  flow 

i ' refers  to  incompressible  flow 
o refers  to  reference  conditions 

T refers  to  turbulent  boundary  layer 

t refers  to  transitional  boundary  layer;  also  refers  to 

stagnation  conditions 

V refers  to  viscous  flow 

w refers  to  the  wall  or  solid  boundary 

x,y  denotes  differention  with  respect  to  x or  y 
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